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A Computer Program Listing Appendix consisting of 95 pages is submitted herewith in hard copy 
and is incorporated herein by reference in its entirety. 

Field of the Invention 



The present invention relates to a techniques for coding and decoding signals used in DSL data 
transmission over wired and wireless systems that use Turbo Codes (TC) and any other receiver soft- 

20 decision extraction technique, such as Low Density Parity Check Codes (LDCP) (these techniques are called 
Turbo-like Codes), with QAM constellations using independent I and Q as modulation for error correction. 
The QAM constellations covered in this application are from 2 QAM up to 16384 QAM. The same 
technique can be use for higher order modulations. In this application we present results using an S-type and 
helicoidal odd-even smile analytical interleaver. An embodiment of the invention pertains particularly to 

25 xDSL systems, as a representative of DMT wired-based systems. The same technique can be used for 

multicarrier wireless DMT systems that using only one transmitter can send multi-channel information. The 
low order constellation insures a minimum acceptable data rate for very low Signal to Noise ratios that is not 
possible with any other known error correcting technique. The preferred mapping used is Gray mapping 
because of its better performance vs. Natural or Ungerboeck mapping, however the techniques disclosed 

30 herein may use Natural or Ungerboeck mapping. 



Background of the Invention 



35 Turbo codes and other receiver soft-decision extraction techniques such as Low Density Parity 

Check Codes, are new and very powerful error control techniques, which allow communication very close 
to the channel capacity. These techniques are called as Turbo-like Codes. Low Density Parity Check Codes 
were introduced in 1962. Turbo codes were introduced in 1993. 



40 A lot of research has been done in the application of Turbo-like Codes in deep space 

communications, mobile satellite/cellular communications, microwave links, paging, in OFDM and CDMA 
architectures. Turbo-like codes outperform all previously known coding schemes regardless of the targeted 
channel The extra coding gain offered by these codes can be used either to save bandwidth or reduce 
power requirements in the link budget Standards based in turbo codes have already been defined or are 

45 currently under investigation. For example: 



• Inmarsat's new multimedia service is based on turbo codes and 16 QAM that allows the user to 
communicate with existing Inmarsat 3 spot beam satellites from a notebook-sized terminal at 
64 kbit/s. 

50 • The Third Generation Partnership Project (3 GPP) proposal for IMT-2000 includes turbo codes 

in the multiplexing and channel coding specification. The IMT-2000 represents the third 
generation mobile radio systems worldwide standard. The 3GPP objective is to harmonize 
similar standards proposals from Europe, Japan, Korea and the United States. 
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• NASA's next-generation deep-space transponder will support turbo codes and implementation 
of turbo decoders in the Deep Space Network is planned by 2003. 

• The new standard of the Consultative Committee for Space Data Systems (CCSDS) in based 
on turbo codes. The new standard outperforms with 1.5 to 2.8 dB the old CCSDS standard 

5 based on concatenated convolutional code and Reed-Solomon code. 

• The new European Digital Video Broadcasting (DVB) standard also adopted turbo codes for 
the return channel over satellite applications. 



xDSL modems are designed to operate between a Central Office CO (or a similar point of 
10 presence) and a customer premises CPE, As such they use existing telephone network wiring between the 
CO and the CPE. There are several modems in this class which function in generally similar manners. 
Typically these modems transmit their signals above the voice band. As such, they are dependent on 
adequate frequency response above voice band. The modem may use a set of transmitters and receivers 
each at different frequencies. Alternatively, a transmitting modem may use the Inverse Discrete Fourier 
15 Transform (IDFT) that sends a signal in different subcarriers at the same time. A corresponding receiver 
uses the Discrete Fourier Transform (DFT) to detect the signal in the different subcarriers. 

The International Telecommunication Union (ITU) has developed Recommendations G.992.1 and 
G.992.2 for the use of ADSL modems in the telephone network. G.992.1 allows the use of Trellis Code 
20 Modulation in the encoder of the transmitter to improve the performance of the ADSL modems. Trellis 
Code Modulation encodes the two least significant information bits using a convolutional encoder that 
provides an extra parity bit for Forward Error Correction. The extra parity bit produces an expansion of 0.5 
bits per dimension in the DMT/QAM modulator. This Forward Error Correction provides an extra coding 
gain of around 3 dB with respect to the un-coded DMT/QAM system. 

25 



Summary of the Invention 

Embodiments of the present invention pertain to the use of Turbo Codes and any other receiver 
30 soft-decision extraction technique such as Low Density Parity Check Codes (these techniques are 

collectively referred to herein as Turbo-like Codes), for DMT systems with QAM modulation in a manner 
that reduces . computing requirements by treating the QAM signal as two AM modulations (one in the I 
dimension and one in the Q dimension) and using the probabilities of the I dimension values as an input (the 
probabilities of the Q dimension are not needed) to the Turbo-like Code. For non-square constellations (8 
35 QAM, 32 QAM, 128 QAM, etc.), we show a method of rectangular constellations with spherical noise; a 
method of creating non-separable I and Q constellations by combining constituent separable I and Q 
constellations, as presented in Provisional Patent Application Serial Number 09/248,099, filed on November 
13, 2000, that also reduces the computational complexity by the order of the square root of the number of 
constellation points (0(N 1/2 )); other methods with lower computational saving such as the diagonal regions 
40 method presented in the Provisional Patent Application Serial Number 09/242,393, filed on October 

20,2000 and Provisional Patent Application Serial Number 09/244,550, filed on October 31, 2000 can also 
be used. 



Further embodiments of the invention improve the Turbo code, and other receiver soft-decision 
45 extraction techniques, by mapping the modulated signal in a manner providing different protection to the 
information bits and the parity bits as a function of the channel used and the performance needed. For 
Additive White Gaussian Noise (AWGN), more protection should be put in the information bits for 
applications with a target BER higher than 10 ~ s , and more protection is preferable put in the parity bits for 
applications with a target BER lower than 10" 8 . More protection to the parity bits than to the information bits 
50 reduces the floor error of the Turbo code below 1 0" 9 . For Impulsive Noise (IN) channels and high 

puncturing, the parity bits should be more protected than the information bits in all cases, obtaining an extra 
protection on the order of 25 dB. 

Further embodiments of the invention employ novel puncturing patterns for each order of a QAM 
55 modulated signal, and a preferred puncturing pattern for each order is identified. The high puncturing used 
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in the Turbo Code makes the Turbo Code and the Low Density Parity Check Code very similar in 
theoretical fundamentals and practical results. 

Embodiments disclosed herein represent the first use of Turbo-like Codes for high order 
5 modulations techniques, up to 16384 QAM. Preferred embodiments achieve a target BER very close to the 
channel capacity for the respective spectral efficiency response. The schemes proposed here, are more 
power efficient than the trellis coded modulation schemes used traditionally with V.32.bis or V,34 
standards. The higher the constellation, the more improvements that this technique provides (i.e. for 16 
QAM the probabilities to compute in a classical turbo code are 16 points, with this technique the 
10 probabilities to compute are 4, so if we use 10 iterations, the computing saving is more than 40 times; for 

128 QAM the computing saving is 10* 1,042/186 = 60 times, using a method of creating non-separable I and 
Q constellations by combining constituent separable I and Q constellations). 

Embodiments disclosed herein apply this technique to DMT xDSL modems as representative 
15 species of a wired-based system, where constellations up to 32,768 QAM are possible, and typically 

constellations of 1024 QAM are used. In the case of the G.992.1, this technique allows the system to work at 
400 kbps with Ej/Nq below 2 dB (assuming 4 QAM and spectral efficiency of 1 bit/s/Hz 4 ksymbols/s and 
100 tones). These embodiments use S-type interleavers and analytical helicoidal odd-even interleavers. The 
attached computer program code listing appendix presents the c code of simulation program, a program for 
20 the S-type interleaver, and a program for the analytical interleaver. 

Brief Description of the Several Views of the Drawings 

Figure 1 shows the System Architecture for variable Symbol rates 

25 

Figure 2 shows the BER curves for the rate 2/4 16 QAM scheme with N=5 12 bits and N-32768 bits S-type 
Figure 3 shows a Coding Scheme 
30 Figure 4 shows a SRC Scheme 

Figure 5 shows a 4 QAM constellation 
Figure 6 shows a 2 AM constellation 

35 

Figure 7 shows the BER for Rate 2/6 4 QAM, N=1024 bits S-type AWGN Channel 
Figure 8 shows the BER for Rate 2/4 4 QAM, N=1024 bits S-type AWGN Channel 
40 Figure 9 shows a 8 QAM constellation 
Figure 10 shows a 4 AM constellation 

Figure 1 1 shows the BER for Rate 2/6 8 QAM, N=1024 bits S-type AWGN Channel 

45 

Figure 12 shows the BER for Rate 4/6 8 QAM, N=1024 bits S-type AWGN Channel 
Figure 13 shows a 16QAM constellation 
50 Figure 14 shows the BER for Rate 1/2 16 QAM, N=272 bits (odd-even) AWGN Channel 
Figure 15 shows the BER for Rate 1/2 16 QAM, N=256 bits S-type AWGN Channel 
Figure 16 shows the BER for Rate 1/2 16 QAM, N=512 bits S-type AWGN Channel 

55 
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Figure 17 shows the BER for Rate 1/2 16 QAM, N=768 bits S-type AWGN Channel 
Figure 18 shows the. BER for Rate 1/2 16 QAM, N-1024 bits S-type AWGN Channel 
5 Figure 19 shows the BER for Rate 3/4 16 QAM, N=4,096 bits S-type AWGN Channel 
Figure 20 shows a 32 QAM constellation 
Figure 21 shows an 8 AM constellation with Gray Mapping 

10 

Figure 22 shows the BER for Rate 3/5 32 QAM, N=6,144 bits S-type AWGN Channel 
Figure 23 shows a 64 QAM constellation 
15 Figure 24 shows a BER for Rate 3/6 64 QAM, N=4,096 bits S-type AWGN Channel 

Figure 25 shows the BER for Rate 4/6 64 QAM, N=4,096 bits S-type AWGN Channel Gray Mapping 
Figure 26 shows an 8 AM constellation with Natural Mapping 

20 

Figure 27 shows the BER for Rate 4/6 64 QAM, N=4,096 bits S-type AWGN Channel Natural Mapping 
Figure 28 shows a 64 QAM constellation: First Partition Level based on Uj = 1 
25 Figure 29 shows a 64 QAM constellation: Second Partition Level based on u x ~ 1, u 2 = 1 

Figure 30 shows a 64 QAM constellation: Third Partition Level based on u t = 1, u 2 = 1, u 3 = 1 
Figure 3 1 shows a 64 QAM constellation: Fourth Partition Level based on u t = 1, u 2 = 1, u 3 = 1, u 4 - 1 

30 

Figure 32 shows a 64 QAM constellation: Fifth Partition Level based on u { = 1, u 2 = 1, u 3 - 1, u 4 = 1, u 5 = 1 
Figure 33 shows a 64 QAM constellation: LLR(u0 
35 Figure 34 shows a 64 QAM constellation: LLR(u 2 ) 
Figure 35 shows a 64 QAM constellation: LLR(u 3 ) 
Figure 36 shows a 64 QAM constellation: LLR(u 4 ) 

40 

Figure 37 shows a 64 QAM constellation: LLR(u 5 ) 

Figure 38 shows a 64 QAM constellation: LLR(u 6 ) 

45 Figure 39 shows the BER for Rate 4/6 64QAM, N=4,096 bits S-type AWGN Channel TCM with four bits 
coded 

Figure 40 shows a 64 QAM constellation with four point subsets 
50 Figure 41 shows a First Partition of the 64 QAM constellation 
Figure 42 shows a Second Partition of the 64 QAM constellation 
Figure 43 shows a Third Partition of the 64 QAM constellation 

55 
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Figure 44 shows a Fourth Partition of the 64 QAM constellation 

Figure 45 shows the BER for Rate 4/6 64 QAM, N=2,048 bits S-type AWGN Channel TCM with two bits 
coded 

5 

Figure 46 shows the BER for Rate 5/6 64 QAM N=5,120 bits AWGN Channel 
Figure 47 shows a 128 QAM constellation 
10 Figure 48 shows a 16 AM constellation with Gray Mapping 

Figure 49 shows the BER for Rate 5/7 128 QAM N=5,120 bits AWGN Channel 
Figure 50 shows a 256 QAM constellation 

15 

Figure 51 shows the BER for Rate 5/8 256 QAM N=5,120 bits AWGN Channel 
Figure 52 shows the BER for Rate 6/8 256 QAM N=6,144 bits AWGN Channel 
20 Figure 53 shows a 512 QAM constellation 

Figure 54 shows the BER for Rate 6/9 5 12 QAM N=6, 144 bits AWGN Channel 
Figure 55 shows a 1024 QAM constellation 

25 

Figure 56 shows the BER for Rate 7/10 1024 QAM N=2,044 bits AWGN Channel 
Figure 57 shows the BER for Rate 12/14 16384 QAM N=3 1,200 bits AWGN Channel 
30 Figure 58 shows a AWGN Channel model 

Figure 59 shows the BER for Rate 2/4 4 QAM N=2,048 bits S-type interleaver AWGN Channel 
Figure 60 shows the BER for Rate 2/4 16 QAM N=2,048 bits S-type interleaver AWGN Channel 

35 

Figure 61 shows the BER for Rate 3/4 16 QAM N=2 ? 046 bits S-type interleaver AWGN Channel 
Figure 62 shows the BER for Rate 4/6 64 QAM N=2,048 bits S-type interleaver AWGN Channel 
40 Figure 63 shows the BER for Rate 5/8 256 QAM N=2,050 bits S-type interleaver AWGN Channel 
Figure 64 shows the BER for Rate 6/8 256 QAM N=2,046 bits S-type interleaver AWGN Channel 
Figure 65 shows the BER for Rate 7/10 1024 QAM N=2 ? 044 bits S-type interleaver AWGN Channel 

45 

Figure 66 shows the BER for Rate 2/4 4 QAM N=2,100 bits Analytical odd-even smile interleaver AWGN 
Channel 

Figure 67 shows the BER for Rate 2/4 16 QAM N=2,100 bits Analytical odd-even smile interleaver AWGN 
50 Channel 

Figure 68 shows the BER for Rate 3/4 16 QAM N=2,100 bits Analytical odd-even smile interleaver AWGN 
Channel 
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Figure 69 shows the BER for Rate 4/6 64 QAM N=2,100 bits Analytical odd-even smile interleaver AWGN 
Channel 

Figure 70 shows the BER for Rate 5/8 256 QAM N=2,100 bits Analytical odd-even smile interleaver 
5 AWGN Channel 

Figure 71 shows the BER for Rate 6/8 256 QAM N=2,100 bits Analytical odd-even smile interleaver 
AWGN Channel 

10 Figure 72 shows the BER for Rate 7/10 1024 QAM N=2,100 bits Analytical odd-even smile interleaver 
AWGN Channel 

Figure 73 shows a 16QAM constellation with indication of xl, x2, x3, x4 and yl, y2, y3, y4 

1 5 Figure 74 shows the performance of the 6/4 rate 64 QAM for the case that the parity bits are more protected 
and for the case that the parity bits are least protected 
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30 



Figure 75 shows a detailed explanation of the encoding and puncturing patterns in a case where all the 
information bits are encoded. 

Figure 76 shows a detail explanation of the encoding and punctured procedures in a case where a subset of 
information bits are encoded. 

Figure 77 shows a process in a transmitter in accordance with embodiments of the invention, and 
Figure 78 shows a process in a receiver in accordance with embodiments of the invention. 

Detailed Description of the Invention 

We herby incorporate by reference the following references, which are representative of the 
conventional knowledge in the field of the invention: 



1 . R. G. Gallager, "Low Density Parity-Check Codes", IRE Trans, on Information Theory, pp. 2 1- 28, 
35 01/62. 

2. C. Berrou, V. Glavieux and P. Thitimajshima, "Near Shannon limit error-correcting coding and 
decoding: turbo-codes", ICC 1993, Geneva, Switzerland, pp, 1064-1070, May 1993. 

3. H. Feldman and D. V. Ramana, "An introduction to Inmarsat's New Mobile Multimedia Service", 
The Sixth International Mobile Satellite Conference, Ottawa, pp. 226-229, June 1999. 

40 4. P. Chaudhury, W. Mohr and S. Onoe, "The 3GPP Proposal for IMT-2000", IEEE Communications 

Magazine, vol 37, no 12, pp. 72- 81. December 1999. 

5. 3GPP Standard "Multiplexing and channel coding; TS 25.212" 

6. C. D. Edwards, C. T. Stelzried, L. J. Deutsch and L. Swanson, "NASAS's deep-Space 
Telecommunications Road Map" TMO Progress Report 42-136, JPL, Pasadena, CA USA. PP. 1-20 

45 February 1999. 

7. R. Pyndiah, A. Picard and A. Glavieux "Performance of Block Turbo Coded 16 QAM and 64 
QAM modulations" Procedings of Globecom 95 pp.1039-1043. 

8. Rauschmayer, Dennis J. "ADSL/VDSL Principles", Macmillan Technical Publishing, 1999. 
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9. ITU G.992.1 "ADSL Transceivers", ITU, 1999. 

10. ITUG.992.2 "Splitterless ADSL Transceivers". ITU 1999. 

11. ITU 1.432 "B-ISDN user-network interface-physical layer specification", ITU , 1993. 

12. Benedetto, Divsalar, Montorsi and F. Pollara, "Serial Concatenation of Interleaved Codes: 

5 Performance Analysis, Design, and Iterative Decoding", The Telecommunications and Data 

Acquisition Progress Report 42-126, Jet Propulsion Laboratory, Pasadena, California, pp. 1-26, 
August 15, 1996. 

13. Benedetto, Divsalar, Montorsi and F. Pollara, "A Soft-Output Maximum A Posteriori (MAP) 
Module to decode parallel and Serial Concatenated Codes", The Telecommunications and Data 

10 Acquisition Progress Report 42-127, Jet Propulsion Laboratory, Pasadena, California, pp. 1-20, 

November 15, 1996. 

14. L. R. Bahl, J. Cocke, F. Jelinek, and J. Raviv, "Optimal Decoding of Linear Codes for Minimizing 
Symbol Error Rate," IEEE Transactions on Information Theory, pp. 284-287, March 1974. 

15. Divsalar and F. Pollara, "Turbo Codes for PCS Applications", Proceedings of ICC'95, Seattle, 
15 Washington, pp. 54-59, June 1995. 

16. D. Divsalar and F. Pollara, "Multiple Turbo Codes", Proceedings of IEEE MILCOM95, San Diego, 
California, November 5-8, 1995. 

17. D. Divsalar and F. Pollara, "Soft-Output Decoding Algorithms in iterative Decoding of Turbo 
Codes," The Telecommunications and Data Acquisition Progress Report 42-124, Jet Propulsion 

20 Laboratory, Pasadena, California, pp. 63-87, February 15, 1995. 

18. Juan Alberto Torres, Frederic Hirzel and Victor Demjanenko, "Forward Error Correcting System 
With Encoders Configured in Parallel and/or Series", International Patent Application Serial No. 
PCT/US99/17369 filed on July 30, 1999. 

25 1. Introduction 

Error correction can be used in communication systems to improve the amount of data that can be 
accurately communicated through a communication channel. Digital Subscriber Line (xDSL) technology is 
an example of communication using error correction. In xDSL, a transmitter and a receiver communicate 

30 over a copper twisted pair using Discrete Multi-Tone (DMT) modulation. In particular, ADSL (Asymmetric 
Digital Subscriber Line) provides multiple carriers in the frequency band between 20 kHz and 1.1 MHz, at 
intervals of 4.1325 kHz. Each carrier is modulated using quadrature amplitude modulation (QAM), in 
which a data value is represented by the phase and magnitude of the carrier. A given phase/magnitude 
combination is referred to as a "symbol", and the set of symbols that can be transmitted on a carrier is 

35 referred to as a "constellation." Typically a constellation contains 2 n symbols; in other words, each symbol 
represents n bits of data. An ADSL channel is theoretically capable of a 1 5 bit constellation (2 15 discrete 
symbols). However, in practice, noise on the line will cause the magnitude and phase of transmitted symbols 
to vary slightly, making it uncertain which symbol a received phase/magnitude combination corresponds to. 
This typically requires the use of a constellation that provides more space between symbols so that 

40 variations in phase and magnitude create less uncertainty. As a result, a given carrier will typically use a 
constellation of less than 15 bits. 
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A type of error correction that is used in xDSL and other communication systems is called Forward 
Error Correction (FEC). Forward error correction involves transmitting one or more error correction bits, 
typically referred to as "parity bits," in conjunction with information bits, so that the information bits 
represented in the signal can be determined solely from the received signal itself using the parity bits, 
5 thereby eliminating the need for retransmission. One simple conventional form of forward error coding 
involves providing a parity bit that represents the least significant bit of a sum of all information bits. Upon 
receiving the bit string, if the least significant bit of the sum of the information bits matches the parity bit, 
the received bits are assumed to be correct. 

10 In a more complex forward error coding system referred to as trellis coding, some input bits are 

coded by a convolutional coder that generates a parity bit for each group of n input bits, and the resulting 
group of n+1 output bits are used by QAM modulation to produce a signal for a digital to analog converter. 
The addition of the parity bit effectively doubles the number of symbols in the constellation, and the value 
of the parity bit effectively places a given symbol in one of two subsets of the constellation, each of which 

1 5 contains every other symbol of the foil constellation. As a result, all symbols of the constellation are 

available for use, but each individual symbol is discriminated within a subset of the constellation in which 
the symbols are twice as far apart as they would be in the full constellation. 

Presently some xDSL modems use a forward error correction method known as Trellis Code 
20 Modulation. Trellis coded modulation (TCM) proposed by Ungerboeck in 1982 is now a well-established 
technique in digital communications. Since its first appearance, TCM has generated a continually growing 
interest, concerning its theoretical foundations as well as its numerous applications, spanning high-rate 
digital transmission over voice circuits, digital microwave radio relay links, and satellite communications. 

25 Turbo-like Codes represent a more recent development in the coding research field (1993), which 

have generated a large interest in the coding community. Turbo Codes are Parallel Concatenated 
Convolutional Codes (PCCC) whose encoder is formed by two (or more) constituent systematic encoders 
joined through one or more interleavers. The input information bits feed the first encoder and, after having 
been scrambled by the interleaver, they enter the second encoder. 

30 

This document describes the performance of Turbo-like Codes used with square and non-square 
high order QAM modulation/demodulation schemes, from 4 QAM up to 1 6384 QAM. For higher order 
modulation the same procedure is applicable. This document also describes the application of this 
modulation schemes to xDSL modems presented in the ITU Recommendations G.992.1 and G.992.2. The 

35 use of squared and non-squared QAM constellations with novel puncturing patterns simplifies the design of 
the equalizer, and for this reason it is possible to utilize squared and non-squared QAM constellations with 
independent I and Q dimensions using blind equalizers. LDPC codes and other receiver soft-decision 
extraction techniques may be employed in a similar manner. The different puncturing patterns allow the 
system to work with low signal-to-rioise signals (in the order of 1 .25 dB) at the expense of reduced data 

40 rates. The Turbo-like Coding and modulation schemes disclosed herein are designed to accommodate a 
variable symbol rate, which could be adapted to a particular channel characteristic. 

The performance of the coding and the modulation schemes disclosed herein are estimated through 
computer simulations, for different interleaver sizes. Bit probabilities are estimated from the received QAM 
45 symbol. Trade offs between power and bandwidth efficiency are alsoaddressed. 



1.1 Objectives 

50 The objectives of using higher order modulations and Turbo-like Codes with independent I and Q are: 

1 . To define a flexible architecture to allow variable symbol rates function of the channel quality. 

2. To estimate the performance of 4 QAM, 8 QAM, 16 QAM, 32 QAM, 64 QAM, 128 QAM, 256 QAM 
512 QAM, 1024 QAM and 16384 QAM schemes for different interleaver sizes and Additive Gaussian 

55 noise (AWGN) and Impulse noise channel (IN). 
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3. To find the mapping technique to limit the constellation expansion to at most 0.5 bit/dimension. 

4. To evaluate the performance for higher order modulation schemes in an AWGN channel and in an 
Impulse Noise Channel. 

5 1.2 System Parameters 

The system parameters are summarized in Table 1. 

Table 1 . System Parameters 

10 



Parameter 


Value 


Target bit error rate (BER) 


10-4, io -5 9 io -6 and 10 "7 


Target Eb/No 


As low as possible 


Interleaver block sizes 


variable to 14 DMT symbols for a latency of 10 ms. 


Modulation 


QAM 


Coding rate 


From 1/3 and up 


Maximum Symbol Rate 


1,024 ksym/s 


Minimum Symbol Rate 


128 ksym/s 


Symbol Rate Step 


4 ksym/s 


Maximum Information Rate 


6,144 kbit/s 


Minimum Information Rate 


160 kbit/s 


Information Rate Step 


32 kbit/s 


Channel 


AWGN and IN 



2. System Architecture. 



2.1 Introduction 

15 

The coding and modulation schemes of the embodiments disclosed herein are designed to 
accommodate variable symbol rates using the same encoder and decoder structure. Therefore, depending on 
the channel quality, the same data rate can be transmitted in a variable bandwidth or a different data rate can 
be transmitted, using the same structure but different puncturing patterns and mapping. 

20 

Figure 1 illustrates functional architectures of a transmitter and a receiver in a communication 
system. At the transmitter, incoming information bits, d, are encoded by a turbo-like coder such as a rate 
1/3 systematic encoder or linear parity bit generator matrix H. The turbo-like encoder outputs the 
information bits d and a parity bit p and parity bit q corresponding to each information bit The puncturing 

25 block performs a puncturing function that determines the actual coding rate of the parity generator matrix H. 
The puncturing function is a simple deletion of some parity bits from the p and q streams which results in 
creation of subsets of the p and q bit streams. No puncturing is applied to the information bits d. A 
mapping function then combines the information bit stream with the subsets of the parity bit streams 
selected by the puncturing pattern and produces a QAM symbol stream by mapping a first subset of the 

30 combined bit streams to an I dimension and mapping a second subset of the combined bit streams to a Q 

dimension. Each puncturing pattern and mapping is identified using a unique word (UW). The puncturing 
pattern and mapping may be varied adaptively. For example, transmission may initially employ the most 
powerful coding scheme that can be used in the available bandwidth. As a function of the Channel State 
Information (with the values of gi and bi for each tone for the case of a G.992.1 or G.992.2 modem), the 

35 puncturing and mapping can be adapted to reduce the transmitted signal bandwidth to a valuethat allows 
communication at the desired BER. The receiver is provided with an identification of the puncturing 
pattern and mapping. 

At the receiver, after synchronization with the transmitter, a particular configuration of the bit 
40 estimator block corresponding to the puncturing pattern and mapping of the transmitter is established. The 
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function of the bit estimator block is to compute the log-likelihood ratios (LLR) for each bit of the received 
QAM symbol in the Turbo-like Code. The outputs from this block are punctured streams of the estimate 
d , p and q bits that are provided to a decoder.. The decoder does not need to know anything about 
the transmitter mapping or modulation scheme. A maximum a posteriori (MAP) algorithm is used in the 
5 decoder to produce decoded bits d. In the case of LDPC the output d is the decoded bit d or an error flag 
indicating that the decoding process has an error. 
2.2 Power vs. Bandwidth. 

The architecture of Figure 1 is flexible enough to select the most bandwidth efficient modulation 
10 given the available signal-to-noise ratio available for a particular channel. An example is shown in Table 2 
for an information data rate of 2,048 kbps. The most powerful coding scheme that could be used is the rate 
2/6 4 QAM. This scheme can achieve a BER of 10" 6 at approximately Ei/Nq = 1,2 dB and requires a symbol 
rate of 3,072 ksym/s. If the transmission occurs over a less noisy channel, them the bandwidth can be 
reduced by using higher spectral efficiency schemes. The minimum symbol rate that can be achieved is 292 
15 ksym/s. 



Table 2. Trade Off Power for Bandwidth for 2,048 kbit/s Information Data Rate 



Spectral efficiency 


Coding Rate 


Symbol Rate 


E b /N 0 




and 






[bits/s/Hz] 


Modulation 


[ksym/s] 


[dB] 


2/3 


2/6 and 4 QAM 


3072 


1.2 


1 


2/4 and 4 QAM 


2048 


2.1 


2 


2/4 and 16 QAM 


1024 


3.9 


3 


3/4 and 16 QAM 


682 


6.0 


4 


4/6 and 64 QAM 


512 


8.8 


5 


5/8 and 256 QAM 


408 


11.5 


6 


6/8 and 256 QAM 


342 


13.8 


7 


7/10 and 1024 QAM 


292 


16.4 



20 The turbo decoder has to operate at a maximum speed of 6, 144 kbit/s. Speeds of 20,048 kbps can 

be achieved with the current FPGA technology 

2.3 Optimization for Specific Applications 

25 Another advantage of this architecture is that the system is re-configurable in order to achieve the 

Quality of Service requested by a particular application. For example, the LDPC code or the turbo code 
modem could be configured to operate in two modes: 

• Mode A: minimum transfer delay for delay sensitive applications. 

30 

In this mode, if a maximum transfer delay is defined, the block for the parity matrix size or the 
interleaver size could be changed in function of the data rate (to an integer number multiple of DMT 
symbols). 

35 • Mode B: maximum coding gain for data transfer applications. 

In this mode, a fixed block parity matrix size or an interleaver size (i.e. 65,536 bits) could be used to 
achieve the best performance. A significant coding gain of almost 2 dB could be achieved at a BER of 
10' 6 if the larger interleaver is used. 

40 

2.4 Why are Turbo-like Codes Different? 

The performance of Turbo-like Codes depends on the delay that is allowed in the turbo-like 
encoder. In the case of Turbo Codes, the elementary decoders rely heavily on the interleaving/de- 
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interleaving process to de-correlate their soft outputs. The larger the interleave^ the better the de- 
correlation that can be achieved, the larger the parity check matrix and the better the performance of the 
LDPC codes. 
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The sensitivity of the Turbo-like Codes to the correlation between data bits can also be noticed in 
higher order modulation schemes when two or more data bits are estimated from the same receiver symbol. 
Different mapping techniques are investigated for a spectral efficiency of four information bits per symbol 
using 64 QAM. Simulation results shows that Gray mapping performs better than Natural mapping, which 
in turn is better than the Ungerboeck set partitioning. A fourth case when two bits where sent un-coded also 
performed badly because of a smaller interleaver size. 

The complexity of the bit estimator block in Figure 1 increases with the constellation size because 
more bits need to be estimated per symbol. The equivalent of this in a trellis coded modulation scheme with 
un-coded bits is that more parallel paths are added to the trellis, therefore there is an increased 
computational effort in the branch metric estimator. 

2.5 Capacity Bound 

The minimum E{/Nq values to achieve the Shannon bound, 4 QAM, 8 QAM, 16 QAM, 32 QAM, 
64 QAM, 128 QAM, 256 QAM, 512 QAM and 1024 QAM bounds for spectral efficiencies form 2/3 up to 7 
bits/s/Hz respectively are as in Table 3 for a BER=10" 5 . 

Table 3 Shannon and QAM bounds. 
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The conversion for E s /No to Ej, /Nq is performed using the following relation 

E b /NofdBJ =E S /N 0 [dB] - 10 log 10 (n) [dB] 
(1) 

where t] is the number of information bits per symbol. 

The required C/Nq given a certain Ef/No can be found using the following relation: 

C/No [dB-Hz] = Eb/No [dB] + 10 logw (Rb) [dB-Hz] 
(2) 

where R}, is the information bit rate. 

For a D-dimension modulation the following formulae are used: 



SNR = 
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SNR = - 



.No r,No 



(4) 



D^- D- 
2 2 

where o^n is the noise variance in each of the D dimension and tj is the number of information bits per 
symbol. From the above relations: 



i-jr ( 2r ^ Eh 
V No 
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3. Coding and Modulation for 2/3 bit/s/Hz Spectral Efficiency 

The only option investigated in this section combines a rate 2/6 coding scheme with 4 QAM. 

3.1 Coding 

A preferred embodiment uses a coding scheme as shown in Figure 3. The encoder is formed by two 
1/2 systematic recursive constituent Concatenated Convolutional Encoders (CCE) that receive two inputs 
every information input bit cycle, adding one parity bit per dimension. We propose the convolution encoder 
format (gl,g2) = (23o,35o) joined through an interleaver and which can be generalized to a N constituent 
systematic encoders with N-l interleavers. 

Figure 4 shows a preferred convolutional encoder, gl and g2 are the values in octal of the forward 
and feedback lines. In Figure 2, the upper side the first line is feedback (value 1), the second line is also a 
feedback (value 2) and the fifth line as well (value 16). The feedback lines sum to a value of 19 in decimal 
or 23 in octal, and is represented as 23o. In Figure 2, the lower side the first line is forward line (value 1), 
the same thing happen with the third (value 4), forth (value 8) and fifth (value 16), this makes the total 
forward line to a value of 29 in decimal or 35 in octal, and is represented as 35o. 

3.2 Puncturing 

In order to obtain a code rate of 2/6, no puncturing is applied. The transmission pattern is given in 



Table 4. 



Table 4. Transmission pattern for Rate 2/6 4 QAM. 



Information bit (d) 


dl 


d2 


parity bit (p) 


PI 


P2 


parity bit (q) 


qi 


q2 


2 AM symbol (I) 


(ui) = (di) 


(ui) = (d 2 ) 


(ui) = (qi) 


2 AM symbol (Q) 


(U2) = (P1 ) 


(U2 ) = (q2 ) 


0*2 ) = (P2 ) 


4 QAM symbol (I,Q) 


(I,Q) = (ui,u 2 ) = (d b pi) 


(W) = (ui,U2) = 

(d 2 ,q 2 ) 


(I,Q) = (ui,U2) = 

(qi>P2) 
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3.3 Modulation 

A 4 QAM scheme is shown in Figure 5. At time k, the symbol u k = (u, k ) is sent through the 
channel and the point A in two dimensional space is received. 

For a 4 QAM constellation with points at -A and A, The E m is: 



■2 A 2 



(6) 



For a rate 2/6 code and 4 QAM, the noise variance is: 
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It is assumed that the time k, u t k modulates the I component and u 2 k modulates the Q component 
for a 4 QAM scheme. 

5 

The symbol u k symbol has the following mapping: u k = (u^, u 2 k ) = (d j , p'); u k+1 = (uf + \ u 2 k+1 ) = 
(d i+ \ q i+1 ); u k+2 = ( Ul k+2 ? u 2 k+2 ) = (q i+ \ p i+2 ). 

Considering two independent Gaussian noises with identical variance o 2 ^ the LLR can be 
10 determined independently for each I and Q. 

At the receiver, the I and Q signals are treated independently in order to take advantage of the 
simpler formulae for the 2 bit-LLR values. 

15 In order to estimate the performance of this scheme, a 2 AM modulation is used, as it is shown in 

Figure 6, instead of 4 QAM modulation. 
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The 4 QAM scheme will achieve a similar performance in terms of bit error rate (BER) at twice the 
spectral efficiency, assuming an ideal demodulator. 

For a rate 2/6 code and a 2 AM scheme as shown in Figure 6, the noise variance is: 

E avI = (l+l)A 2 /2 = A 2 (8) 
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That is the same that the 4 QAM of equation (7) 



30 3.4 Bit Probabilities 

For an A WGN channel the following expressions need to be evaluated: 



35 



LLR( u ))=\og * 
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The above LLRs are used as inputs to the turbo encoder. There is no need to compute the 4 LLRs 
for all symbols because I and Q signals are treated independently. Also the simulations for the 2 bit-LLR 
40 values are reduced to one term each. 
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3.5 Simulations Results 

Figure 7 shows the performance of a turbo code which 1024 bit interleaves using S-type 
interleaver. 

The target BER of 10 " 7 for a 1,024 information bit interleaver can be achieve at Eft/No = 1.25 dB 

4. Coding and Modulation for 1 bit/s/Hz Spectral Efficiency 

Two options are investigated in this section. The first scheme combines a rate 2/4 coding scheme 
with 4 QAM. The second scheme investigated combines 2/6 coding scheme with 8 QAM. 

4. 1 Option 1 : Rate 2/4 Turbo Code and 4 QAM 

4.1.1 Coding 

The coding scheme is shown in Figure 3. 

The two systematic recursive codes (SRC) used are identical and are defined in Figure 4, The code 
is described by the generating polynomials 35o and 23o. 

4.1.2 Puncturing 

In order to obtain a code rate of 2/4, every other bit of the parity bits p and q from Figure 3 are 
punctured. The puncturing pattern is given in Table 5. 

Table 5. Puncturing and Mapping for Rate 2/4 4 QAM. 



Information bit (d) 


dl 


d 2 


parity bit (p) 


PI 




parity bit (q) 




q2 


2 AM symbol (I) 


(ui) = (di) 


(ui) = (d 2 ) 


2 AM symbol (Q) 


(U2) = (P1 ) 


0*2 ) = to ) 


4 QAM symbol (I,Q) 


(I,Q) = (ui f u 2 )-(di,pi) 


(I,Q) = (ui,u 2 ) = (d2,q2) 
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4.1.3 Modulation 
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A 4 QAM scheme is shown in Figure 5. At time k, the symbol u k - (u^, u 2 k ) is sent through the 
channel and the point in two dimensional space is received. 



For a 4 QAM constellation with points at -A and A, The is: 

4(A 2 + A 2 )_, 



Eav 



- = 2 A 2 



For a rate 2/4 code and 4 QAM, the noise variance is: 
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It is assumed that the time k, k modulates the I component and u 2 k modulates the Q component 
for a 4 QAM scheme. 
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The symbol u k symbol has the following mapping: u k = (u, k , u 2 k ) = (d 1 , p'); u k+1 = (u,^ 1 , u 2 k+I ) = 
(d' +, ,q' +I ). 

Considering two independent Gaussian noises with identical variance c 2 ^ the LLR can be 
determined independently for each I and Q, 

At the receiver, the I and Q signals are treated independently in order to take advantage of the 
simpler formulae for the 2 bit-LLR values. 

In order to estimate the performance of this scheme, rate 1/2 turbo codes a 2 AM modulation is 
used, as it is shown in Figure 6, instead of 4 QAM modulation. 

The 4 QAM scheme will achieve a similar performance in terms of bit error rate (BER) at twice the 
spectral efficiency, assuming an ideal demodulator. 



For a rate 1/2 code and a 2 AM scheme as shown in Figure 6, the noise variance is: 
E avI = (l+l)A 2 /2 =A 2 
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That is the same that the 4 QAM of equation (13) 
4.1.4 Bit Probabilities. 

For an AWGN channel the following expressions need to be evaluated: 

f i / . ^ f 
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The above LLRs are used as inputs to the turbo encoder. There is no need to compute the 4 LLRs 
for all symbols because I and Q signals are treated independently. Also the simulations for the 2 bit-LLR 
values are reduced to one term each. Due to the puncturing, with one in two parity bits being transmitted, 
the expected performance will be lower when compared with the non-punctured scheme. 

For an Impulse Noise channel, for a determined period of time (the duration of the impulse) the 
E^/Nq ratio changes to a worse value, and the errors occurs. The Reference Level of Ei/Nq is set to the 
value of Ef/NQ needed by the un-coded QAM constellation to reach the value of BER of 10" 7 . 

4.1.5 Simulations Results 
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Figure 8 shows the performance of a turbo code which 1024 bit interleaves using S-type 
interleaves 
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The target BER of 10 " 7 for a 1 ,024 information bit interleaver can be achieve at Ef/No = 2. 1 dB 
4.2 Option 1: Rate 1/3 Turbo Code and 8 QAM 
5 4.2.1 Coding 

The coding scheme is shown in Figure 3. 

The two systematic recursive codes (SRC) used are identical and are defined in Figure 4. The code 
10 is described by the generating polynomials 35o and 23o. 

4.2,2 Puncturing 

No puncturing is applying in this case. The pattern is given in Table 6, 

15 

Table 6. Pattern and Mapping for Rate 1/3 8 QAM. 



Information bit (d) 


dl 


parity bit (p) 


PI 


parity bit (q) 


qi 


4 AM symbol (I) 


(ui,u 2 ) = (di, pi) 


2 AM symbol (Q) 


(«3)=(qi) 


8 QAM symbol (I , Q) 


(I,Q) = (dl,Pl,qi) 
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4.2,3 Modulation 

The 8 QAM scheme used is shown in Figure 9. At time k, the symbol u k = (u^, u 2 k , u 3 k ) is sent 
through the channel and the point in two dimensional space is received. 



For an 8 QAM constellation with points at -3A, -A, A and 3 A, in the I dimension and -A, and A in 
25 the Q dimension. 
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The Em, _Q ls: 



Em_Q ~ A 
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so: 



Because the noise has to be spherical, the total value of E^ will be the addition of these two values, 
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Eav =Eav_I + E<av J) = 6 A 2 
(20) 

The value of the noise variance has to be the same in both dimensions and its value is: 
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It is assumed that the time k, ^ k and u 2 k modulates the I component and u 3 k modulates the Q 
component for an 8 QAM scheme. 
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The symbol u k has the following mapping: u k = (u^, u 2 k , u 3 k ) = (d 1 , p', q'). The parity bits are 
mapped to the least protected bits of the QAM symbol. Note that k denotes the symbol time index and i the 
information bit time index. 

Considering two independent Gaussian noises with identical variance o 2 ^ the LLR can be 
determined independently for each I and Q. 

At the receiver, the I and Q signals are treated independently in order to take advantage of the 
simpler formulae for the LLR values. The mapping of the information bit is made to the most protected bit 
in each dimension (uj k for the I signal and u 3 k for the Q signal in Figure 9). 

In order to estimate the performance of this scheme, rate 1/3 turbo codes and 4 AM modulation (as 
it is shown in Figure 10) is used in the I dimension, and a 2 AM modulation is used for the Q dimension, (as 
it is shown in Figure 6), instead of 8 QAM modulation. 

4.2.4 Bit Probabilities 

For an AWGN channel the following expressions need to be evaluated: 
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The above LLRs are used as inputs to the turbo encoder. There is no need to compute the 8 LLRs 
for all symbols because I and Q signals are treated independently. Also the simulations for the LLR values 
are reduced to two terms for the 4 AM modulation and one term for the 2 AM modulation. Due to the 
puncturing, the expected performance will be lower when compared with the non-punctured scheme. 

4.2,5 Simulations Results 

Figure 1 1 shows the performance of a turbo code, which used 1024 bit interleavers using S-type 
interleaver. 

The target BER of 10 " 7 can be achieve at Ejy/No = 3.1dB. 
5. Coding and Modulation for 2 bit/s/Hz Spectral Efficiency 
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Two options are investigated in this section. The first scheme combines a rate 4/6 coding scheme 
with 8 QAM. The second scheme investigated combines 1/2 coding scheme with 16 QAM. These schemes 
have the advantage of a very efficient implementation without any significant compromise in performance. 

5.1 Option 1: Rate 4/6 Turbo Code and 8 QAM 

5.1.1 Coding 

The coding scheme is shown in Figure 3. 

The two systematic recursive codes (SRC) used are identical and are defined in Figure 4. The code 
is described by the generating polynomials 35o and 23o. 

5.1.2 Puncturing 

In accordance with an embodiment of the invention, the puncturing pattern given in Table 7 is used 
in order to obtain a rate 4/6 code. 
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Table 7. Puncturing and Mapping for Rate 4/6 8 QAM. 



Information bit (d) 




d 2 


d3 


d 4 


parity bit (p) 


PI 








parity bit (q) 
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4 AM symbol (I) 


(ui,U2) = (di, pi) 


0n,U2) = (d3, q3) 


2 AM symbol (Q) 


(u 3 ) = (d2) 


(u 3 ) = (d 4 ) 


8 QAM symbol (I,Q) 


(I,Q) = (dl,Pl,d2) 


(I,Q) = (d3>q3><i4) 



5.1.3 Modulation 
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The 8 QAM scheme used is shown in Figure 9. At time k, the symbol u k = (u^, u 2 k ? u 3 k ) is sent 
through the channel and the point r k in two dimensional space is received. 



For an 8 QAM constellation with points at -3A, -A, A and 3 A, in the I dimension and -A, and A in 
the Q dimension. 
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The Egy j is: 
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Because the noise has to be spherical, the total value of E m will be the addition of these two values, 
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£av _ F'CivJ + Eav_Q - 6 A " 
(27) 



The value of the noise variance has to be the same in both dimensions and its value is: 
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It is assumed that the time k, Uj k and u 2 k modulates the I component and u 3 k modulates the Q 
component for an 8 QAM scheme. 
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The symbol u k has the following mapping: u k = (u ( k , u 2 k , u 3 k ) = (d 1 , p l , d 1+l ) and for the consecutive 
symbol: u k+1 = (u 1 k+1 ,u 2 k+1 ) u 3 k+1 ) = (d 1+2 , q 1+2 > d i+3 ). The parity bits are mapped to the least protected bits of the 
QAM symbol. Note that k denotes the symbol time index and i the information bit time index. 

Considering two independent Gaussian noises with identical variance o 2 ^ the LLR can be 
determined independently for each I and Q. 

At the receiver, the I and Q signals are treated independently in order to take advantage of the 
simpler formulae for the LLR values. The mapping of the information bit is made to the most protected bit 
in each dimension (u^ for the I signal and u 3 k for the Q signal in Figure 9). 

In order to estimate the performance of this scheme, rate 4/6 turbo codes and 4 AM modulation (as 
it is shown in Figure 10), is used in the I dimension, and a 2 AM modulation is used for the Q dimension, 
(as it is shown in Figure 6), instead of 8 QAM modulation. 

5.1.4 Bit Probabilities 

For an A WGN channel the following expressions need to be evaluated: 
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The above LLRs are used as inputs to the turbo encoder. There is no need to compute the 8 LLRs 
for all symbols because I and Q signals are treated independently. Also the simulations for the LLR values 
are reduced to two terms for the 4 AM modulation and one term for the 2 AM modulation. Due to the 
puncturing, the expected performance will be lower when compared with the non-punctured scheme. 

5.1.5 Simulations Results 

Figure 12 shows the performance of a turbo code, which used 1024 bit interleavers using S-type 
interleaver. 

The target BER of 10 " 7 can be achieve at the following values of Eft/No ~ 5.5 dB 
5.2 Option 2: Rate 1/2 Turbo Code and 16 QAM 
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5.2.1 Coding 

The coding scheme is shown in Figure 3. 

The two systematic recursive codes (SRC) used are identical and are defined in Figure 4, The code 
is described by the generating polynomials 35o and 23o. 

5.2.2 Puncturing 

In order to obtain a rate 1/2 code, every other bit of the parity bits p and q from Figure 3 are 
punctured. The puncturing pattern is given in Table 8. This is the traditional puncturing pattern using 
alternating parity bits in each cycle. 

Table 8. Puncturing and Mapping for Rate 1/216 QAM. 



Information bit (d) 




d 2 


parity bit (p) 


PI 




parity bit (q) 






4 AM symbol (I) 


(ui,u2) = (di ,pi ) 


4 AM symbol (Q) 


(u 3 ,U4) = (d2,q2) 


16 QAM symbol (I > Q) 


( I , Q) = (ui, U2 , U3, U4 ) = (di , pi, d 2 , q2 ) 
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5.23 Modulation 

A 16 QAM scheme is shown in Figure 13. At time k, the symbol u k = (Uj k , u 2 k , u 3 k , u 4 k ) is sent 
through the channel and the point A in two dimensional space is received. 



For a 16 QAM constellation with points at -3A, -A,A and 3 A, The Eay is: 
r _ 4(A 2 + A 2 )+8( A 2 + 9A 2 )+4(9 A 2 + 9A } ) _ in , 2 

Em 77 JVA 
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For a rate 1/2 code and 16 QAM, the noise variance in each dimension is: 
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It is assumed that the time k ? u t k and u 2 k modulates the I component and u 3 k and u 4 k modulates the 
Q component for a 16 QAM scheme. 

The symbol u k symbol has the following mapping: u k = (u^, u 2 \ u 3 \ u^ = (d\ p 1 , d 1+1 ,q The 
parity bits are mapped to the least protected bits of the QAM symbol. Note that k denotes the symbol time 
index and i the information bit time index. This means a puncturing of one in two parity bits. 

Considering two independent Gaussian noises with identical variance a 2 ^ the LLR can be 
determined independently for each I and Q. 

At the receiver, the I and Q signals are treated independently in order to take advantage of the 
simpler formulae for the 4 bit-LLR values. The mapping of the information bit is made to the most 
protected bit in each dimension (uj k for the I signal and u 3 k for the Q signal in Figure 13). 

In order to estimate the performance of this scheme, rate 1/2 turbo code and 4 AM modulation is 
used (as it is shown in Figure 10), instead of 16 QAM modulation. 
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The 16 QAM scheme will achieve a similar performance in terms of bit error rate (BER) at twice 
the spectral efficiency, assuming an ideal demodulator. 

For a rate 1/2 code and a 4 AM scheme as shown in Figure 6, the noise variance is: 

E av i= (I+9)A 2 /2 =5A 2 
(34) 
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5.2.4 Bit Probabilities 
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For an AWGN channel the following expressions need to be evaluated: 
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25 The above LLRs are used as inputs to the turbo encoder. There is no need to compute the 16LLRs 

for all symbols because I and Q signals are treated independently. Also the simulations for the 4 bit-LLR 
values are reduced to 2 terms each. Due to the puncturing, with one in two parity bits being transmitted, the 
expected performance will be lower when compared with the non punctured scheme. 

30 5.2.5 Simulations Results 

Figure 14 shows the performance of a turbo code that uses the odd-even block helical inter leaver of 
272 bits. Figure 1 5 shows the performance of a turbo code that uses the S-type interleaver for 256 bits. 



35 The S-type interleaver performs better than the odd-even block helical interleaver with 

approximately 1 dB at a BER=10" 6 in this case. Therefore, the rest of the simulations used the S-type 
interleaver. 
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Figures 16 to 18 show the simulation results for 512, 768 and 1024 bit interleaves using S-type 
interleaver. 



5 A search for the best S-type interleaves was performed. For each interleaver size, the maximum 

S-factor was found and four different interleaves were generated using a random seed. Four simulations 
were generated for a target BER of 10" 5 . The best result decided the preferred interleaver, complete 
simulations were started for BERs down to 10 " 7 . 

10 In the case of N=768 } a further interleaver search is required in order to avoid the flattening of the 

BER curve which started to show below 10" 6 in Figure 17. 

5.2.6 Recommended Solution 

15 The rate 1/2 16 QAM scheme described in this chapter achieves the target BER at less than 1.5 dB 

from capacity. The implementation in hardware is feasible and it can be used at very high data rates. 

The target BER of 10" 7 can be achieve at the following values of El/No = 6.8 dB for N=256 
information bits and Eb/Ng= is 4.5 dB for N=1024 information bits. 

20 

6. Coding and Modulation for 3 bit/s/Hz Spectral Efficiency 

Three options are investigated in this section. The first scheme combines a rate 3/4 coding scheme 
with 16 QAM. The second scheme combines a rate 3/5 coding scheme with 32 QAM and the third scheme 
25 combines a rate 3/6 coding scheme with 64 QAM. 

6.1 Option 1 Rate 3/4 Turbo Code and 16 QAM 

6. LI Coding 

30 

The coding used in this scheme is shown in Figure 3. The two systematic recursive codes (SRC) 
used are identical and is defined in Figure 4. The code is described by the generating polynomials 35o and 
23o. 

35 6.1.2 Puncturing 

In accordance with an embodiment of the invention, in order to obtain a rate 3/4 code, the 
puncturing pattern used is given in Table 9. 

40 Table 9. Puncturing and Mapping for Rate 3/4 1 6 QAM 



Information bit (d) 




d 2 


d 3 


d 4 


d 5 


d6 


parity bit (p) 




P2 










parity bit (q) 














4 AM symbol (I) 


(dl,d 2 ) 


(d 4 , d 5 ) 


4 AM symbol (Q) 


(<>3 > P2 ) 


(d6,q5) 


16 QAM symbol (I , Q) 


(I,Q)" (di > d 2 ,d 3 ,p2) 


(I,Q)- (d 4 ,d5,d 6 ,q 5 ) 



6.1.3 Modulation 



45 A 16 QAM scheme is shown in Figure 13. It is assumed that at time k u t k and u 2 k modulates the I 

component and u 3 k and u 4 k modulates the Q component of a 16 QAM scheme. 

In order to estimate the performance of this scheme, rate 3/4 turbo codes and 4 AM modulation are 

used. 
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For a rate 3/4 code and 16 QAM, the noise variance in each dimension is: 
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The puncturing and mapping scheme is shown in Table 9 for 6 consecutive information bits that 
are encoded into 8 coded bits, therefore two 16 QAM symbols. 



For the 4 AM scheme as shown in Figure 13, the noise variance is: 



A2 f2xh5x Eb ] J _lO A2 (E^ J 
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(41) 



k No J 
6.1,4 Bit Probabilities 

For each received symbol, the bit probabilities are computed as described in equations (36) to (39). 
6.L5 Simulation Results 

Figure 19 shows the simulation results for 6,144 information bits with S-type interleaves 
A BER of 10" 7 can be achieved after 8 iterations at Ef/Ng = 5 55 <*B f° r N = 6,144 information bits. 
6.2 Option 2 Rate 3/5 Turbo code and 32 QAM 
25 6.2,1 Coding 

The coding used in this scheme is shown in Figure 3. The two systematic recursive codes (SRC) 
used are identical and are defined in Figure 4. The code is described by the generating polynomials 35o and 
23o. 



6.2,2 Puncturing 

In accordance with an embodiment of the invention, in order to obtain a rate 3/5 code, the 
puncturing pattern used is given in Table 10. 

Table 10. Puncturing and Mapping for Rate 3/5 32 QAM 



Information bit (d) 


dl 


d 2 


d3 


d 4 


d 5 


d6 


parity bit (p) 


PI 






P4 






parity bit (q) 














8 AM symbol (I) 


(dl, d3 , pi) 


(d4> d& qs) 


4 AM symbol (Q) 


(<*2 > 32 ) 


(<*5 > P4 ) 


32 QAM symbol (I,Q) 


(I,Q)= (di,d3,pi,d 2 , q2) 


( I , Q) = (d^dfc q5 ,d 5 , p 4 ) 
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6.2.3 Modulation 

The asymmetrical 32 QAM modulation scheme considered is shown in Figure 20. At time k, the 
symbol u k = (u^, u 2 \ u 3 k , u 4 \ u 5 k ), is sent through the channel and the point i* in two dimensional space is 
received. 
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It is assumed that at time k u,\ u 2 k and u 3 k modulate the I component and u 4 k and u 5 k modulate the Q 
component of a 32 QAM scheme. 

In the I dimension, two information bits and one parity bit are transmitted using Gray mapping and 
8 AM as shown in Figure 2 1 . 

The E m i is: 

fl+9 + 25 + 49) A\ 



Eav i ~ 



2\A 2 



The Eav gis: 



_ (l + 9)A ! _ 

Em:_Q " ->A 



(42) 



(43) 



so: 



Because the noise has to be spherical, the total value of will be the addition of these two values, 

Eav =Eav_l +E av _Q = 26A 2 
(44) 

The value of the noise variance has to be the same in both dimensions and its value is: 

<jn = Eav | | = 26 A 2 



2r]E b ) „ A 2x3x£ b f_26 J E^'* 
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(45) 



The puncturing and mapping scheme is shown in Table 10 for 6 consecutive information bits that 
are encoded into 10 coded bits, therefore two 32 QAM symbols, 

6.2,4 Bit Probabilities 

For an AWGN channel, the following expressions need to be evaluated for the I dimension: 
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LLR( U " 3 ) = \og 



j=; ^ (Tat 



= log 
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For the Q dimension, only LLRJu^ and LLR(u 5 k ) are required to be evaluated: 



(48) 
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6.2.5 Simulation Results 

Figure 22 shows the simulation results for 6,144 information bits with S-type interleaves 

A BER of 10** can be achieved after 8 iterations at Ei/Nq = 6.5 dB for N = 6,144 information bits. 

6.3 Option 2 Rate 3/6 Turbo code and 64 QAM 

6.3.1 Coding 

The coding used in this scheme is shown in Figure 3. The two systematic recursive codes (SRC) 
used are identical and are defined in Figure 4. The code is described by the generating polynomials 35o and 
23o. 



25 6.3,2 Puncturing 

In order to obtain a rate 3/6 code, the puncturing pattern used is given in Table 1 1 . This is the 
traditional puncturing pattern using alternating parity bits in each cycle. 



Table 11. Puncturing and Mapping for Rate 3/6 64 QAM 



Information bit (d) 


<*1 


<*2 


d3 


d 4 


d 5 


d6 


parity bit (p) 


PI 




P3 




P5 




parity bit (q) 




q2 




q4 




% 


8 AM symbol (I) 


(dl,d2,Pl) 


(d4, d5, q4) 


8 AM symbol (Q) 


(d3, P3 , <\2) 


(<*6> P5 > %) 


64 QAM symbol (I , Q) 


(I,QKdl,d2>Pb d 3> P3^2) 


(I ,Q)= (d 4 ,d 5 , q4,d6, P5,q6 ) 



6.33 Modulation 
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Let us consider 64 QAM modulation scheme as it is shown in Figure 23. At time k, the symbol u k 
= (Uj k , u 2 k , u 3 k , u 4 k , u 5 k , u 6 k ), is sent through the channel and the point I* 1 " in two dimensional space is received. 

It is assumed that at time k u 3 k , u 2 k and u 3 k modulates the I component and u 4 k , u 5 k and u 6 k modulates 
the Q component of a 64 QAM scheme. 

For 64 QAM constellations with points at -7A, -5A, -3 A, -A, A, 3 A, 5A, 7A The E av is: 

E av = (8(49+25+9+l)+8(25+49+49+9+49+l)+8(25+9+25+l)+8(9+l)) A 2 / 64= 42 A 2 (51) 

For a rate 3/6 code and 64 QAM, the noise variance in each dimension is: 
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In order to estimate the performance of this scheme, when rate 3/6 turbo codes and 8 AM 
1 5 modulation is used, as shown in Figure 2 1 . 

The 64 QAM scheme will achieve a similar performance in terms of bit error rate (BER) at twice 
the spectral efficiency, assuming an ideal modulator. 



For the 8 AM scheme shown in Figure 21, the noise variance is 



ctn~ E[\ a 2 J 



27] E b )\-2lA 2 ^ 2xl ' SxEb 



V No J V No 

This is identical to the variance for the 64 QAM scheme. 



(53) 



The puncturing and mapping scheme is shown in Table 1 1 for 6 consecutive information bits that 
are encoded into 12 coded bits, therefore two 64 QAM symbols. 

6.3.4 Bit Probabilities 

For an AWGN channel, the following expressions need to be evaluated for the I dimension: 
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An identical computation effort is required for the Q dimension, the I k being replaced with the Q k 
demodulated value in order to evaluate LLR(u 4 k ), LLR(u 5 k ) and LLR(u 6 k ). 

6.3.5 Simulation Results 

Figure 24 shows the simulation results for 6,144 information bits with S-type interleaves A BER of 
10" 7 can be achieved after 8 iterations at Ej/Nq = 6.1 dB for N = 4,096 information bits. This result is 0.3 dB 
worse than the performance of the rate 3/4 16 QAM scheme. 

6.4 Recommended Solution 

Comparing Figures 19, 22 and 24, the recommended solution for 3 bits/Hz spectral efficiency is the 
rate 3/4 turbo code 16 QAM scheme. 

This scheme achieves a target BER of 10* 7 at E^/Nq =5.8 dB for N=6,144 information bits. 

7. Coding and Modulation for 4 bit/s/Hz Spectral Efficiency 

This section investigated four schemes for transmission of 4 information bits in a 64 QAM symbol. 
The mapping used for 64 QAM constellation has a very significant impact on the performance of these 
schemes. The first scheme uses independent I and Q mapping and also uses Gray mapping in each 
dimension. The second scheme uses independent I and Q mapping, but natural mapping in each dimension. 
The third scheme used a conventional trellis coded modulation approach based on Ungerboeck set 
partitioning of the 64 QAM set. This partitioning technique splits the constellation in sub-constellations 
with increased Euclidian distance. In these schemes all the information bits are coded. The fourth scheme 
used the same conventional trellis coded modulation approach based on Ungerboeck set partitioning of the 
64 QAM set. However, only two information bits are encoded by a rate 1/2 code. The four coded bits select 
a sub-partition of four points. The other two information bits, which are sent uncoded, identify the 
transmitted point. 

7.1 Option 1 - Rate 4/6 64 QAM with independent I and Q and with Gray Mapping 

7.1.1 Coding 

The coding scheme is shown in Figure 3. The two systematic recursive codes (SRC) used are 
identical and are defined in Figure 4. The code is described by the generating polynomials 35o and 23o. 

7.1.2 Puncturing 

In accordance with an embodiment of the invention, in order to obtain a rate 4/6 code, the 
puncturing pattern used is shown in Table 12. 

Table 12. Puncturing and Mapping for Rate 4/6 64 QAM Option 1 
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Information bit (d) 


dl 




d3 


d 4 


parity bit (p) 


PI 








parity bit (q) 






43 




8 AM symbol (I) 


(dl,d2,pi) 


8 AM symbol (Q) 


0*3* d 4> Q3) 


64 QAM symbol (I , Q) 


(I,QHdl s d2 9 pi,d3,d4,q3) 
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7.1.3 Modulation 



A 64 QAM scheme is shown in Figure 23. Gray mapping was used in each dimension as shown in 



Figure 21. 



Four information bits are required to be sent using a 64 QAM constellation. 
For a rate 4/6 code and 64 QAM, the noise variance in each dimension is 



ai = E m \^ \ =42 A 2 

No 



2x4xEb 
No 



V 



= 5.25 A — 
\N 0 



(57) 



Gray mapping was used in each dimension as shown in Figure 2 1 . 

15 The puncturing and mapping scheme is shown in Table 12 for 4 consecutive information bits that 

are encoded into 6 coded bits, therefore one 64 QAM symbol. 

The turbo encoder with the puncturing presented in Table 12 is a rate 4/6 turbo code which in 
conjunction with 64 QAM gives a spectral efficiency of 4 bits/s/Hz. 

20 

Considering two independent Gaussian noises with identical variance o 2 ^ the LLR can be 
determined independently for each I and Q. 

It is assumed that at time k u^, u 2 k and u 3 k modulates the I component and u 4 k , u 5 k and u 6 k modulates 
25 the Q component of the 64 QAM scheme. 

At the receiver, the I and Q signals are treated independently in order to take advantage of the 
simpler formulae for the LLR values. 

30 7.1.4 Bit Probabilities 

From each received symbol, the bit probabilities are computed as described in equations (46) (47) 
and (48) for I dimension. An identical computation effort is required for the Q dimension, the I k being 
replaced with the Q k demodulated value in order to evaluate LLR^*), LLR(u 5 k ) and LLR(u 6 k ). 



7.1.5 Simulation Results 

Figure 25 shows the simulation results for 4,096 information bits with S-type interleaver. 

A search for the best S-type interleaver was performed. For each interleaver size, the maximum S- 
factor was found and four different interleavers were generated using a random seed. Four simulations were 
generated for a target BER of 10" 5 , The best result decided the preferred interleaver. Using the preferred 
interleaver, complete simulations were started for BERs down to 10" 7 . 

A BER of 10" 7 can be achieved after 8 iterations at E^No = 8.3 dB. 

7.2 Option 2 - Rate 4/6 64 QAM with independent I and Q and Natural Mapping 
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7.2.1 Coding 

The proposed coding scheme is shown in Figure 3. The two systematic recursive codes (SRC) 
5 used are identical and are defined in Figure 4. The code is described by the generating polynomials 35o and 
23o. 



7.2.2 Puncturing 

10 In accordance with an embodiment of the invention, in order to obtain a rate 4/6 code, the 

puncturing pattern used is shown in Table 13. 



Table 13. Puncturing and Mapping for Rate 4/6 64 QAM Option 2 



Information bit (d) 


dl 


d2 


d 3 


d 4 


parity bit (p) 


PI 








parity bit (q) 






q3 




8 AM symbol (I) 


(dl,d2,pi) 


8 AM symbol (Q) 


( d 3> d 4> q3) 


64 QAM symbol (I,Q) 


(I,QMdi ? d 2 ,pi,d3>cl4,q3) 



7.2.3 Modulation 



A 64 QAM scheme is shown in Figure 23. Natural mapping was used in each dimension as shown 
in Figure 26. Four information bits are required to be sent using a 64 QAM constellation. This is equivalent 
20 to a rate 2/3 coding scheme. 

For a rate 4/6 code and 64 QAM, the noise variance in each dimension is 



25 



30 



2 _ v ] 2rjE b 

No , 



No J 



= 5.25 A 2 



\No 



E b i (58) 



Gray mapping is used in each dimension as shown in Figure 21. 

The puncturing and mapping scheme is shown in Table 13 for 4 consecutive information bits that 
are encoded into 6 coded bits, therefore one 64 QAM symbol. 

The turbo encoder with the puncturing presented in Table 13 is a rate 4/6 turbo code which in 
conjunction with 64 QAM gives a spectral efficiency of 4 bits/s/Hz. 



Considering two independent Gaussian noises with identical variance cr^, the LLR can be 
35 determined independently for each I and Q. 

It is assumed that at time k u^, u 2 k and u 3 k modulates the I component and u 4 k , u 5 k and u 6 k modulates 
the Q component of the 64 QAM scheme. 

40 At the receiver, the I and Q signals are treated independently in order to take advantage of the 

simpler formulae for the LLR values. 

7.2.4 Bit Probabilities 



45 From each received symbol, the bit probabilities are computed as described in equations (46) (47) 

and (48) for I dimension. An identical computation effort is required for the Q dimension, the I k being 
replaced with the Q k demodulated value in order to evaluate LLR(u 4 k ), LLR(u 5 k ) and LLR(u 6 k ). 
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7.2.5 Simulation Results 
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Figure 27 shows the simulation results for 4,096 information bits with S-type interleaver. 

A search for the best S-type interleaver was performed. For each interleaver size, the maximum S- 
factor was found and four different interlevers were generated using a random seed. Four simulations were 
generated for a target BER of 10" 5 . The best result decided the preferred interleaver. Using the preferred 
interleaver, complete simulations were started for BERs down to 10~ 7 . 

A BER of 10" 7 can be achieved after 8 iterations at Ej/No = 10.5 dB. 

7.3 Option 3 - Trellis Coded Modulation with 4 bits Coded 

7.3.1 Coding 

In this scheme, all four information bits are coded by a rate 4/6 code. Only two parity bits are 
transmitted. The six bits produced select a point in the 64 QAM constellation. The coding scheme is shown 
in Figure 3. The two systematic recursive codes (SRC) used are identical and are defined in Figure 4. The 
code is described by the generating polynomials 35o and 23o. 

7.3.2 Puncturing 

In accordance with an embodiment of the invention, in order to obtain a rate 2/3 code, the 
puncturing pattern used is shown in Table 14. 

Table 14. Puncturing and Mapping for Rate 4/6 64 QAM Option 3 



Information bit (d) 




d2 


d 3 


d 4 


parity bit (p) 


PI 








parity bit (q) 
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64 QAM symbol (I , Q) 


(I,QMdl,d 2 , d 3 , d 4 ,pi ? q 3 ) 



30 7.3.3 Modulation 

A trellis coded modulation scheme is employed. The six bits divide the 64 QAM constellation 
based on increased Euclidean distance. 

35 The 64QM constellation is partitioned in subsets as shown in Figures 28 to 32. 

The puncturing and mapping scheme is shown in Table 14 for 4 consecutive information bits that 
are encoded into 6 coded bits, therefore one 64 QAM symbol. 

40 7.3.4 Bit Probabilities 

For an AWGN channel, the following expressions need to be evaluated for each received symbol 
before the turbo decoding process can start. 
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LLR( U \) = log 
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(51) 



The II R* - P, II represents the squared Euclidian distance between the received point R* at the time 
k and a collation poinfp, The actual set in which u, is 0 or 1 are shown m F lg ures 33-38. 

15 7.3.5 Simulation Results 

Tlis scheme required much higher computational effort than previous options and would be 

A BER of 10" 7 can be achieved after 8 iterations at El/No = 115 dB - 
n a n r tinn 4 - Trellis Pvted Modulati ™ with 2 hits Coded. 
7.4.1 Coding 

Tn thisscheme only two information bits are coded by a rate 2/4 code. The other twc .information 
35oand23o. 
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7,4.2 Puncturing 



In accordance with an embodiment of the invention, in order to obtain a rate 1/2 code, every other 
parity bit is punctured as shown in Table 15. 

Table 15. Puncturing and Mapping for Rate 1/2 64 QAM Option 4 



Information bit (d) 




d2 


d 3 


d 4 


parity bit (p) 


PI 








parity bit (q) 










64 QAM symbol (I , Q) 


(I,Q) = (di,pi,d 3 , d2,q2>d 4 ) 



10 



7.4.3 Modulation 

A trellis coded modulation scheme is employed. The four coded bits divide the 64 QAM 
constellation based on increased Euclidean distance. 



The 64QM constellation shown in Figure 39 is partitioned by the four coded bits in 16 subsets with 
15 four point each. The two uncoded bits will select one of the four points of the subset. Figures 41 to 44 show 
the four steps of the partitioning process. Each 16 points constellation subset can now be further partitioned 
as show in Figures 4 1 to 44. 

The puncturing and mapping scheme is shown in Table 15 for 4 consecutive information bits that 
20 are encoded into 6 coded bits, therefore one 64 QAM symbol. 

7.4.4 Bit Probabilities 

For an AWGN channel, the following expressions need to be evaluated for each received symbol 
25 before the turbo decoding process can start. 
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Each summation in equations (64) (65) (66) and (67) is over 32 points 

The || R k - P t || represents the squared Euclidian distance between the received point R k at the time 
k and a constellation point The actual set in which Uj is 0 or 1 are shown in Figures 41-44. 

5 

7.4.5 Simulation Results 

Figure 45 shows the simulation results for 4,096 bit blocks with 2,048 bit S-type interleaves 

10 A BER of 10~ 7 can be achieved after 8 iterations at E}/Nq = 11.5 dB. 

7.5 Recommended Solution 

Comparing the complexity of each scheme and their performance, the recommended solution for 4 
15 bit/s/Hz spectral efficiency is Option 1, the rate 4/6 turbo code 64 QAM scheme independent I and Q with 
gray mapping in each dimension. This achieves the target BER at Ej/Nq = 8.3 dB. For N=4096 information 
bits. 

8. Coding and Modulation for 5 bit/s/Hz Spectral Efficiency. 

20 

This section investigated three schemes all use independent I and Q modulation. The first scheme 
combines a rate 5/6 coding scheme with 64 QAM. The second scheme combines a rate 5/7 coding scheme 
with 128 QAM. The third scheme combines a rate 5/8 coding scheme with 256 QAM. 

25 8.1 Option 1 - Rate 5/6 Turbo Code and 64 QAM. 

8.1.1 Coding 

The coding scheme is shown in Figure 3. The two systematic recursive codes (SRC) used are 
30 identical and are defined in Figure 4. The co de is described by the generating polynomials 35o and 23o. 

8.1.2 Puncturing 

In accordance with an embodiment of the invention, in order to obtain a rate 5/6 code, the 
35 puncturing pattern used is shown in Table 16. 



Table 16. Puncturing and Mapping for Rate 5/6 64 QAM Option 1 



Information bit (d) 




d2 


d3 


d 4 


d5 


d6 


d 7 


dg 


d 9 


dio 


parity bit (p) 


PI 




















parity bit (q) 






















8 AM symbol (I) 


(dl,d2,pi) 


(d& d 7 > q6) 


8 AM symbol (Q) 


(d3, d4, ds) 


(dg, d9, djo) 


64 QAM symbol (I , Q) 


<I,QMdl,d2,pi,d3,d4,d5) 


( I ,Q)=(ds,d7> q6>d8, d^dio) 



40 8.1.3 Modulation 

A 64 QAM scheme is shown in Figure 23 . 

Five information bits are required to be sent using a 64 QAM constellation. This is equivalent to a 
45 rate 5/6 coding scheme. 

For a rate 5/6 code, the noise variance in each dimension is 
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The puncturing and mapping scheme is shown in Table 16 for 10 consecutive information bits that 
are encoded into 12 coded bits, therefore two 64 QAM symbols. The turbo encoder with the puncturing 
presented in Table 16 is a rate 5/6 turbo code which in conjunction with 64 QAM gives a spectral efficiency 
of5bits/s/Hz. 

Considering two independent Gaussian noises with identical variance a 2 ^ the LLR can be 
determined independently for each I and Q. 

It is assumed that at time k Uj k , u 2 k and u 3 k modulates the I component and u 4 k , u 5 k and u 6 k modulates 
the Q component of the 64 QAM scheme. 

At the receiver, the I and Q signals are treated independently in order to take advantage of the 
simpler formulae for the LLR values. 

8. 1.4 Bit Probabilities 

From each received symbol, the bit probabilities are computed as described in equations (54) (55) 
and (56) for I dimension. An identical computation effort is required for the Q dimension, the I k being 
replaced with the Q k demodulated value in order to evaluate LLR(u 4 k ), LLR(u 5 k ) and LLR(u 6 k ). 

8. 1 .5 Simulation Results 

Figure 46 shows the simulation results for 5,120 information bits (1,204 QAM symbols) with S- 
type interleaves 

The high puncturing rate makes the iterative decoding process to converge very slowly, showing a 
flattening of the BER curve. Therefore this option is not considered acceptable. 

8.2 Option 2 - Rate 5/7 Turbo Code and 128 QAM. 

8.2.1 Coding 

The coding scheme is shown in Figure 3. The two systematic recursive codes (SRC) used are 
identical and are defined in Figure 4. The co de is described by the generating polynomials 35o and 23 o. 

8.2.2 Puncturing 



40 In accordance with an embodiment of the invention, in order to obtain a rate 5/7 code, the 

puncturing pattern used is shown in Table 17. 

Table 17. Puncturing and Mapping for Rate 5/7 128 QAM Option 2 



Information bit (d) 


dl 


<*2 


d3 


d 4 


d5 


parity bit (p) 


PI 










parity bit (q) 
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16 AM symbol (I) 


(d 3 , d4, ds s q3) 


8 AM symbol (Q) 


(dl. d 2, Pi) 


128 QAM symbol (I , Q) 


(I,QHd3»d4 s <*5> Q3> di,d 2? pi) 



45 



8.2.3 Modulation 
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Let us consider the 128 QAM scheme as shown in Figure 47. It is assumed that at time k the 
symbol u k = (u^, u 2 k , u 3 k , u 4 k , u 5 k , u 6 k ) is sent though the channel It is assumed that at time k the symbol U] k , 
u 2 k , u 3 k and u 4 k modulate the I component and u 5 \ u 6 k and u 7 k modulate the Q component of a 128 QAM 
scheme. 

In the I dimension, three information bits and one parity bits are transmitted using Gray mapping 
and 16 QAM as show in Figure 48. 

The average energy The E av I for the 16 QAM shown in Figure 48 is: 

The EqvJ is: 

(T+9 + 25 + 49 + 81 + 121 + 169 + 225;^* . 2 
E^_j = " 85/| 



(69) 



The Egy gis: 

_ (T+9 + 25+49)^ „ 2 

Eav_Q ~ Z[ A 



(70) 



so: 



Because the noise has to be spherical, the total value of E^ will be the addition of these two values, 

Eav =EavJ+EavJ2= 106 A 2 
(71) 

The value of the noise variance has to be the same in both dimensions and its value is: 
( 2t]E h 



2 . 

On ' 



No 



106 a 



2 \ 2x5x Eb 
No 



= 10.6 A 2 



Eb 
\Noj 



(72) 
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The puncturing and mapping scheme shown in Table 17 is for 5 consecutive information bits that are coded 
into 7 coded bits, therefore, one 128 QAM symbol The turbo encoder is a rate 5/7 turbo code, which in 
conjunction with 128 QAM, gives a spectral efficiency of 5 bits/s/Hz. 

8,2.4 Bit Probabilities 

The 16 QAM symbol is defined as u k = (Ui k , u 2 k , u 3 k , u 4 k ), where u^ is the most significant bit and u 4 k is 
the least significant bit . The following set can be defined. 





1. 


bit-l-is-0 = 


{ A 0 ,A X , A 2 , A 3 , A 4 , A 5 , A 6 , A 7 } 
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2, 


bit-l-is-1 = 


{ A 8 ,A 9 , A 10 , A n , A 12 , A ]3 , A 14 , A 15 } 




3. 


bit-2-is-0 -- 


{ A 0 >A! , A 2 , A 3 , A 8 , A 9 , A 10 , A n } 




4. 


bit-2-is-l = 


{ A 4 ,A 5 , A 6 , A 7 , A 12 , A 13 , A 14 , A 15 } 




5. 


bit-3-is-0 = 


{ A 0 ,A^ , A 4 , A 5 , A s , A 9 , A 12 , A 13 } 




6. 


bit-3-is-l = 


{ A 2 ,A 3 , A 6 , A 7 , A 10 , A u , A 14 , A 15 } 
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7, 


bit-4-is-0 = 


{ A 0 ,A 2 , A 4 , A 6 , A 8 , A 10 , A 12 , A H } 




8. 


bit-4-is-l = 


{ A! ,A 3 , A 5 , A 7 , A 9 , A u , A 13 , A 15 } 



From each received symbol, R , the bit probabilities are computed as follows: 
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f ( j 

I ex P -ttI**-^ 



LLR( u l 1 ) = log 



*> 2 1 



I exp(-Alk-^: 



(73) 
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LLR( u \)=\og 



LLR( u k sMog 



I txpi-j-rlit-A, 



(74) 



(75) 



LLRU) = \o% 



I exp . — l^-^J 



yAjebit-4-i s-0 



(76) 



Similar formulae apply for the 8 AM symbol defined as u k = (u 5 k , u 6 k , u 7 k ), where u 5 k is the most 
significant bit and u 7 k is the least significant bit. 

10 8.2.5 Simulation Results 

Figure 49 shows the simulation results for 5,120 information bits (1,204 QAM symbols) with S- 
type interleaves 

1 5 The high puncturing rate makes the iterative decoding process to converge very slowly, showing a 

flattening of the BER curve. 

8.3 Option 3 - Rate 5/8 Turbo Code and 256 QAM 

20 8.3.1 Coding 

The coding scheme is shown in Figure 3. The two systematic recursive codes (SRC) used are 
identical and are defined in Figure 4. The co de is described by the generating polynomials 35o and 23o. 

25 8.3.2 Puncturing 

In accordance with an embodiment of the invention, in order to obtain a rate 5/8 code, the 
puncturing pattern used is shown in Table 18. 
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Table 18. Puncturing and Mapping for Rate 5/8 256 QAM Option 3 



Information bit (d) 


dl 


d 2 


d 3 


d 4 


d5 


d6 


d 7 


d8 


d 9 


dlO 


parity bit (p) 


PI 








P5 






P8 






parity bit (q) 






Q3 






q6 








qio 


16 AM symbol (I) 


(di, d2,d3,pi) 


(d& d7 , dg, q6) 


16 AM symbol (Q) 


(d4> ds, q3 ,ps) 


(d9,di0, P8 ,qio) 


256 QAM symbol (I , Q) 


(I,QHdi 42 43 >P1 4445 A3 , 
P5) 


(I,Q)=(d6 47 48 >q6 49410 >P8 > 
Q10) 



8.3.3 Modulation 

35 Let us consider the 256 QAM scheme as shown in Figure 50 
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For a 256 QAM constellation with points at -15 A, -13 A, -1 1 A, -9A, -7A, -5A, -3A, -A, A, 3A, 5A, 
7A, 9A, 11 A, 13A, 15A.E av is: 

E av = (8(225+ 1 69+ 121+8 1+49+25+9+1 )+8(225+ 1 69+225+ 1 2 1+225+8 1+225+49 
5 +225+25+225+9+225+ 1 )+8( 1 69+ 121+1 69+8 1+1 69+49+ 1 69+25+1 69+9+ 1 69+ 1 ) (77) 

+8(121+81+121+49+121+25+121+9+121+l)+8(81+49+81+25+81+9+81+l) 
+8(49+25+49+9+49+l)+8(25+9+25+l)+8(9+l)) A^ / 256 = 170 A 2 

It is assumed that at time k the symbol u k = (u t k , u 2 k , u 3 k , u 4 k , u 5 k , u 6 k , u 7 k , u 8 k ) is sent though the channel. 
10 It is assumed that at time k the symbol u 5 k , u 2 k , u 3 k and u 4 k modulate the I component and u 5 k , u 6 k , u 7 k and 
u s k modulate the Q component of a 256 QAM scheme. 

For a rate 5/8 code and 256 QAM, the noise variance is: 



15 al- 



20 



25 



A No J { No 



= 17 A 



J* (78) 

No) 



In order to study the performance of this scheme, a rate 5/8 turbo code and a 16 AM is used as 
describe in Figure 50. The 256 QAM scheme will achieve a similar performance in terms of bit error rate 
(BER) at twice the spectral efficiency, assuming an ideal demodulator. 

The E av for the 16 QAM shown in Figure 48 is: 

E av i = (1+9+25+49+81+121+169+225) A 2 / 8 = 85 A 2 
(79) 

The noise variance is: 



CTN 



l - R J & Y =85 £*H£ & Y .„J Ml]' (80) 

No J {No 



TO 



30 This is identical to the variance for 256 QAM scheme. 

The pimcturing and mapping scheme shown in Table 18 is forlO consecutive information bits that 
are coded into 16 encoded bits, therefore, one 256 QAM symbol. 

35 The turbo encoder is a rate 5/8 turbo code, which in conjunction with 256 QAM, gives a spectral 

efficiency of 5 bits/s/Hz. 

8.3.4 Bit Probabilities 

40 The 16 AM symbol is defined as u*^ (Uj k , u 2 \ u 3 k , u^, where Uj k is the most significant bit and u 4 k 

is the least significant bit . The following set can be defined. 





1. 


bit-l-is-0 




{ A 0 ,A| , A 2 , A 3 , A 4 , A 5 , A 6 , A 7 } 




2. 


bit-l-is-1 




{ A 8 ,Ag , A 10 , A n , A i2 , A 13 , A 14 , A l5 } 
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3. 


bit-2-is-0 




{ A 0 y A Y , A 2 , A 3 , A 8 , A 9 , A i0 , A n } 




4. 


bit-2-is-l 




{ A 4 ,A 5 , A 6 , A 7 , A 12 » A 13 , A 14 , A 15 } 




5. 


bit-3-is-0 




{ A 0 )A X , A 4 , A 5 , A 8 , A 9 , A 12 ? A 13 } 




6. 


bit-3-is-l 




{ A 2 ,A 3 j A 6 , A 7 , A 10 , A u , A 14 , A 15 } 




7. 


bit-4-is-O 




{ A 0 ,A 2 , A 4 , A 6 , A 8 , A l0 , A 12 , A 14 } 
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8. 


bit-4-is-l 




{ Aj ,A 3 , A 5 , A 7 , A 9 , A n , A 13 , A 15 } 



From each received symbol, R k , the bit probabilities are computed as follows: 



37 



LLRtf,) = tog 



A,<=bit-I-is-l \ *<JN J 



LLR(u\)=\o% 



I ex P Yj-iU'-a, 



A^bU-2-is-} 



'<<TN 



A^btt-Us-l V ^^jV 



I 

yAj^bit-3-is-O 



ex Pi-^-r«-R 
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LLRtfj= log 
8.3.5 Simulation Results 



A t tzb i t-4-i s-1 



(81) 



(82) 



(83) 



(84) 



Figure 51 shows the simulation results for 5,120 information bits (1,204 QAM symbols) with S- 
10 type interleaves 

A search for the best S-type interleaver was performed. For each interleaver size, the maximum S- 
factor was found and four different interleaves were generated using a random seed. Four simulations were 
generated for a target BER of 10* 5 . The best result decided the preferred interleaver. Using the preferred 
15 interleaver, complete simulations were started for BER down to 10* 7 . A BER of 10~ 7 can be achieved after 8 
iterations at Ef/No = 11-8 dB. 

8.4 Recommended Solution 

20 The recommended solution for 5 bit/s/Hz spectral efficiency is the rate 5/8 turbo code 256 QAM 

scheme independent I and Q with gray mapping in each dimension. This achieves the target BER at E^/Nq 
= 1 1.8 dB. For N = 5 ? 120 information bits. 



9. Coding and Modulation for 6 bit/Hz Spectral Efficiency 

This section investigated two schemes all use independent I and Q modulation. The first scheme 
combines a rate 6/8 coding scheme with 256 QAM. The second scheme combines a rate 6/9 coding scheme 
with 512 QAM. 

30 9.1 Option 1 - Rate 6/8 Turbo Code and 256 QAM 

9.1.1 Coding 

The coding scheme is shown in Figure 3. The two systematic recursive codes (SRC) used are 
35 identical and are defined in Figure 4. The code is described by the generating polynomials 35o and 23o. 



9. 1 .2 Puncturing 
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In accordance with an embodiment of the invention, in order to obtain a rate 6/8 code, the 
puncturing pattern used is shown in Table 19. 



Table 19. Puncturing and Mapping for Rate 6/8 256 QAM Option 1 

5 



Information bit (d) 






d3 


d 4 


d5 


d6 


parity bit (p) 


PI 












parity bit (q) 














16 AM symbol (I) 


(Ah d 2 , d3, Pi) 


16 AM symbol (Q) 


(d 4 , d5 , d& q4) 


256 QAM symbol (I, Q) 


( 1 9 Q ) = (di,d2, d3, pi, d4,d5, d 6 , q 4 ) 



9.1.3 Modulation 
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Let us consider a 256 QAM scheme as shown in Figure 50. It is assumed that at time k the symbol 
u k = (u^, u 2 k , u 3 k , u 4 k , u 5 k , u 6 k , u 7 k , u 8 k ) is sent though the channel. It is assumed that at time k the symbol u^, 
u 2 k , u 3 k and u 4 k modulate the I component and u 5 k , u 6 k , u 7 k and u s k modulate the Q component of a 256 QAM 
scheme. 



For a rate 6/8 code and 256 QAM, the noise variance is: 



No ) 



No 



No) 



(85) 



The E av in the I and Q dimension is the 16 QAM shown in Figure 47 is: 

E av i = (1+9+25+49+81+121+169+225) A 2 / 8 - 85 A 2 
(86) 

The noise variance is: 
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cn Ea 



2flR 

No 



= S5A 



2{ 2x3x£b 



= 14.16^1^ 



No J \No 
This is identical to the variance for 256 QAM scheme. 



(87) 



The puncturing and mapping scheme shown in Table 19 is for 6 consecutive information bits that 
are coded into 8 coded bits, therefore, one 256 QAM symbol. The turbo encoder is a rate 6/8 turbo code, 
which in conjunction with 256 QAM ? gives a spectral efficiency of 6 bits/s/Hz. 

9.1.4 Bit Probabilities 

The 16 AM symbol is defined as u k = (u^ , u 2 k , u 3 k , u 4 k ), where u^ is the most significant bit and 
u 4 k is the least significant bit . The following set can be defined. 





1. 


bit-l-is-0 = 


{ ?A] , A 2 , A 3 , A 4 , A 5 , A 6 , A 7 } 




2. 


bit-l-is-1 = 


{ A 8 ,A 9 , A 10 , A H , A 12 > A 13 , A I4 , A 15 } 
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3. 


bit-2-is-0 = 


{ A 0 9^1 ? A 2 , A 3 , A g , A 9 , A 10 , A u } 




4. 


bit-2-is-l = 


{ A 4 ,A 5 , A 6 , A 7 , A 12 , A 13 , A 14 , A 15 } 




5. 


bit-3-is-0 = 


{ A 0 ,Aj , A 4 , A 5 , A 8 , A 9 , Aj 2 , A ]3 } 




6. 


bit-3-is-l = 


{ A 2 ,A 3 , A 6 , A 7 , A 10 , A n , A 14 , A 15 } 




7. 


bit-4-is-0 = 


{ A 0 ,A 2 , A 4 , A 6 , A s , A i0 , A ]2 , A 14 } 
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8. 


bit-4-is-l = 


{ Aj ,A 3 , A 5 , A 7 , A 9 , A u , A 13 , A 15 } 
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From each received symbol, R k , the bit probabilities are computed as follows: 



LLR( u k ,)=\og 



'A; 



LLR( u k 2 )=\og 



f 1 ^ 
I exp -_]|^-^|| 



(88) 



A r <Ebu-2-i s-1 



£<Jn 



(89) 



LLR( u k s) = \og 



Z ex p(-y^lk-^ 



r r / v 



(90) 



LLR( u k 4 ) = log 



Aj^bit-4-i s-0 



(91) 
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9.L5 Simulation Results 



Figure 52 shows the simulation results for 6 ? 144 information bits (1,204 QAM symbols) with S- 
15 type interleaves 

A search for the best S-type interleaver was performed. For each interleaver size, the maximum S- 
factor was found and four different interleaves were generated using a random seed. Four simulations were 
generated for a target BER of 10~ 5 . The best result decided the preferred interleaver. Using the preferred 
20 interleaver, complete simulations were started for BER down to 10" 7 . 

A BER of 10 -7 can be achieved after 8 iterations at EfJNo = 14.2 dB, 

9.2 Option 2 - Rate 6/9 Turbo Code and 5 12 QAM 

9.2.1 Coding 
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The coding scheme is shown in Figure 3. The two systematic recursive codes (SRC) used are 
identical and are defined in Figure 4. The code is described by the generating polynomials 35o and 23o. 

9.2.2 Puncturing 

In accordance with an embodiment of the invention, in order to obtain a rate 6/9 code, the 
puncturing pattern used is shown in Table 20, 

Table 20. Puncturing and Mapping for Rate 6/9 512 QAM Option 2 



Information bit (d) 


d] 


d2 


c*3 


d 4 


d5 


d6 


d 7 




d 9 


dlO 


dll 


dl2 


parity bit (p) 


PI 








P5 








P9 
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parity bit (q) 


q 3 - 


q 7 - qn | - 


32 AM symbol (I) 


(d 1; d2,d3,d4 ) pi) 


(dj, dg , d% dio , q7) 


16 AM symbol (Q) 


(<*5> d 6> P5 A3) 


(dn, di2,qn ? P9) 


512 QAM symbol (I , Q) 


ftQM<*l ^2 43 > PI A d 6 ,P5 > Q3) 


(I,Q)=(d7 , dg , dg ,dio ,q7 ,dn,di2 ,qi 1 

>P9) 



15 



20 



9,2.3 Modulation 

Let us consider a 5 12 QAM scheme as shown in Figure 52. It is assumed that at time k the symbol 
5 u k = (u L k , u 2 k , u 3 k , u 4 k , u 5 \ u 6 k , u 7 k , u 8 k , u 9 k ) is sent though the channel. It is assumed that at time k the symbol 
Uj k , u 2 k , u 3 k , u 4 k and u 5 k modulates the I component and u 6 k , u 7 k , u s k and U9 k modulates the Q component of a 
512 QAM scheme. 

In the I dimension, four information bits and one parity bit are transmitted using Gray mapping and 
10 32 AM as shown in Figure 53 . 

The 32 constellation points have amplitudes: -31A, -29A ? -27A, -25A, -23A, -21 A, -19A, -17A, - 
ISA, -13A, -1 1A, -9A, -7A, -5A, -3A, -A, A, 3A, 5A, 7A, 9A, 1 1A, 13A, 15A, 17A, 19A, 21A, 23A, 25A, 
27 A, 29A,31A.TheE avl is: 



The EqwJ is: 

EavJ = 341 A 

The £ m; gis: 

(T+9 + 25 + 49+81+121 + 169 + 225)^ 2 _ _ 2 
Eav Q $5 A 



(92) 



(93) 



Because the noise has to be spherical, the total value of E^ will be the addition of these two values, 



so: 
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Eav =EavJ+EavJ2 = 426 A 2 

The value of the noise variance has to be the same in both dimensions and its value is: 



CTjV I~ E a 



No J V No 



341 
12 ' 



Ml 

No 



(94) 



(95) 



The puncturing and mapping scheme shown in Table 20 is for 12 consecutive information bits that 
are coded into 1 8 coded bits, therefore, one 512 QAM symbol The turbo encoder is a rate 6/9 turbo code, 
30 which in conjunction with 512 QAM, gives a spectral efficiency of 6 bits/s/Hz. 

9.2.4 Bit Probabilities 

From each received symbol, R k , the bit probabilities are computed as follows: 



LLR( u )) = \og 



2 «xp\-^TT\lt-A 



LLR( U \) = log 



E ex pf-yrlk-^l 

AitbU-2-is-l V ^tFJV 



(96) 



(97) 
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LLR(u)) = \ou 



LLR( u " 4 )=log 



{ ( I W 

I exp - — l^-^ll 



2 ex P \-Tjl*-4 



2ai' 



A,e btt-4-Js-J V Z (?y J_ 



LLRful ) = \og 



2 J-^\\r>-a,\\) 



X ex p 

s-0 



(98) 



(99) 



(99) 



9.2.5 Simulation Results 
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Figure 54 shows the simulation results for 6,144 information bits (1,204 QAM symbols) with S- 
type interleaves 

A BER of 1 0" 6 can be achieved after 8 iterations at EtfNo = 1 3 .5 dB. A BER of 1 0~ 7 can be 
achieved after 8 iterations at EyNo = 14.5 dB. 

9.3 Recommended Solution 

The recommended solution for 6 bit/s/Hz spectral efficiency is the rate 6/8 turbo code 256 QAM 
scheme. This scheme achieves a BER of 10~ 7 at Ef/NQ = 14.2 dB. 

10, Coding and Modulation for 7 bit/Hz Spectral Efficiency 

This section investigated one scheme that use independent I and Q modulation. The scheme 
combines a rate 7/10 coding scheme with 1024 QAM. 

10.1 Coding 

The coding scheme is shown in Figure 3. The two systematic recursive codes (SRC) used are 
identical and are defined in Figure 4. The co de is described by the generating polynomials 35o and 23o. 

10.2 Puncturing 

In accordance with an embodiment of the invention, in order to obtain a rate 7/10 code, the 
puncturing pattern used is shown in Table 2 1 . 

Table 21. Puncturing and Mapping for Rate 7/10 1024 QAM 



Information bit (d) 




d2 


d 3 


d 4 


d5 


d6 


dy 


dg 


d 9 


dlO 


dll 


dl 


dl 


dl 


























2 


3 


4 


parity bit (p) 


PI 










P6 










Pll 








parity bit (q) 


























qi 

3 




32 AM symbol (I) 


(dl,d2,d3, pi ,q3) 


(d8,d9,dio, dn . q8 ) 


32 AM symbol (Q) 


(d4,d5,d6, d7,p6) 


(di2,di3,di4,pn ,qi3) 
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1024 QAM symbol (I,Q) 



(I,Q) = (di,d 2 , d 3s pi, q 3 ,d 4 ,d 5 , d6, d 7 , 
P6) 



(IQ) = 

(d8>d9 3 dl04ll a q8 g dl2»dl3»di4»Pll>qi3) 
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10.3 Modulation 

Let us consider a 1024 QAM scheme as shown in Figure 55. It is assumed that at time k the symbol 
u k = (u^, u 2 \ u 3 \ u 4 k , u 5 \ u 6 k , u 7 k , u 8 k u 9 k , u^*) is sent though the channel It is assumed that at time k the 
symbol u^, u 2 k , u 3 k , u 4 k and u 5 k modulates the I component and u 6 k , u 7 \ u 8 \ u 9 k and u 10 k modulates the Q 
component of a 1024 QAM scheme. 

For a 1024 QAM constellation with points at -3 1A, -29A, -27A, -25A, -23A, -21A, -19A, -17A, - 
ISA, -13A, -1 1A, -9A, -7A, -5A, -3A, -A, A, 3A, 5A, 7A, 9A, 11A, 13A, ISA, 17A, 19A, 21A, 23A, 25A, 
27A, 29A,31A,E av is: 

E n = 8(961+841+729+625+529+441+361+289+225+1 69+ 12 1 +8 1+49+25+9+ 1 )+8(96 1+841+96 1 

+729+961+625+961+529+961+441+961+361+961+289+961+225+961+169+961+121+961+81+961 

+49+961+25+961+9+961+l)+8(841+729+841+625+841+529+841+441+841+361+841+289+841+ 

225+841+169+841+121+841+81+841+49+841+25+841+9+841+l)+8(729+625+729+529+729+441 

+729+361+729+289+729+225+729+169+729+121+729+81+729+49+729+25+729+9+729+l)+8(625+ 

529+625+441+625+361+625+289+625+225+625+169+625+121+625+81+625+49+625+25+625 (100) 

+9+625+l)+8(529+441+529+361+529+289+529+225+529+169+529+121+529+81+529+49+529+25+ 

529+9+529+ 1)+8(441+361+44 1+289+441+225+441+169+441+121+441+81+441+49+441+25+441 

+9+441+l)+8(361+289+361+225+361+169+361+121+361+81+361+49+361+25+361+9+361+l)+ 

8(289+225+289+169+289+121+289+81+28^-49+289+25+289+9+289+l)+8(225+169+225+121 

+225+81+225+49+225+25+225+9+225+l)+8(169+121+169+81+169+49+169+25+169+9+169+l) 

+8(121+81+121+49+121+25+121+9+121+l)+8(81+49+81+25+81+9+81+l)+8(49+25+49+9+49+l) 

+8(25+9+25+ 1 )+8(9+ 1 )) A 2 /1024 = 682 A 2 



For a rate 7/10 code and 512 QAM, the noise variance is: 



No ) 



No 



= 4S.7A 2 



\No) 



(101) 



The 1024 QAM constellation is a product of two 32 constellation with points at: -3 1 A, -29A, -27 A, 
-25A, -23A, -21 A, -19A, -17A, -ISA, -13A, -11A, -9A, -7A, -5A, -3 A, -A, A, 3A, 5A, 7 A, 9A, 1 1A, 13A, 
ISA, 17A, 19A, 21 A, 23 A, 25A, 27A, 29A, 31 A. The 32 QAM E avJ is: 

E av = (1+9+25+49+81+121+169+225+289+361+441+529+625+729+841+961) A 2 / 16 = 341 A 2 
(102) 



For a rate 7/10 code and 1024 QAM, the noise variance is: 



CTN ~ Ea 



= 48.7,4' 



\N 0 



(103) 



v \ 

V No ) \ No 

This is identical to the variance for 1024 QAM scheme. 

The puncturing and mapping scheme shown in Table 21 is for 14 consecutive information bits that 
are coded into 20 coded bits, therefore, two 1024 QAM symbols. The turbo encoder is a rate 7/10 turbo 
code, which in conjunction with 1024 QAM, gives a spectral efficiency of 7 bits/s/Hz. 

10.4 Bit Probabilities 

The 32 AM symbol is defined as u k = (u^, u 2 \ u 3 k , u 4 \ u 5 k ), where Uj k is the most significant bit and 
u 5 k is the least significant bit. The following set can be defined. 

1. bit-l-is-0={ Ao, Ai , A2 , A3, A4, A5 , A6, A7, Ag, A9,Aio, An , Aj2 , A\$ , A14 , A15 } 

43 



10 



2. bit-l-is-1 = { Aj6 ,A\1 , A]8 , A19 , A20 , A21 , A22 , A23 , A24 ,A25 , A26 > A27 , A28 , A29 , A30 , A31 } 

3. bit-2-is-0 = { Aq , Ai ? A2 , A3 , A4 , A5 , A$ , A7 , A\$ ,A\7 , Aig , A19 , A20 , A21 , A22 > A23 } 

4. bit-2-is-l - { As 9 A9 , A10 , An, A12 , A13 , A14 , A15, A24 ,A25 , A26 , A27 , A28 , A29 , A30 , A31 } 

5. bit-3-is-0 = { Ao , Ai , A2 , A3, Ag a A9 , A10 , An, Ai6 All , M% > A i9, A24 ,A25 , A26 , A27 } 

6. bit-3-is-l = { A4 , A5 , A6 , A7, Ai2 , A13 , A\4 , A\ 5 , A20 , A21 , A22 , A23, A28 , A29 , A30 , A3 1 } 

7. bit-4-is-0 = { Aq , A] , A4 , A5 , Ag , A9 , An , A13 , Ai6 ,Aj7 , A20 , A21 , A24 , A25 , A28 , A29 } 

8. bit-4-is-l *= { A2 , A3 , A6 , A7 , A10 , A] 1 , A14 , A15 , Ais ,Ai9 , A22 , A23 , A26 , A27 , A30 , A3 1 } 

9. bit-5-is-0 = { Ao , A2 , A4 , A6 , Ag , A10 , A12 , A14 , A\$ ? A]g , A20 , A22 , A24 , A26 , A28 , A30 } 

10. bit-5-is-l = { Ai , A3 , A5 , A7 , A9 , An , A13 , A15 , A17 ,Aj9 , A21 , A23 , A25 , A27 , A29 , A31 } 

From each received symbol, R k , the bit probabilities are computed as follows: 
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LLR(u) ; = log 



R k -Ai 



LLR( U \) = log 



A, ebii-2-is-J 



1 



H exp 

\^Aj£bit-2-* s-0 



1 



(104) 



(105) 
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LLR(u) ; = log 



2 exp 

A.ebifS-is-J 



•<JN 



Ajebtt-3-is-O \ ^GN 



R k -Aj\ 



n 



LLR( u l ,)=\og 



LLRtf,) = log 



2 «pf-rrl« i -A|)' 

A^bil-4-is-l V )_ 

[jebit-4-is-O V ) 

4,eiit-S-is-l \ J 



^ Aj^bit-5-i s-Q 



\*-As\ 



(106) 



(107) 



(108) 
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10.5 Simulation Results 



Figure 56 shows the simulation results for 2 ? 044 information bits with S-type interleaver. 



A search for the best S-type interleaver was performed. For each interleaver size, the maximum S- 
factor was found and four different interleavers were generated using a random seed. Four simulations were 
30 generated for a target BER of 1 0' 5 . The best result decided the preferred interleaver. Using the preferred 
interleaver, complete simulations were started for BER down to 10" 7 . 

A BER of 10" 7 can be achieved after 8 iterations at Ej/Nq = 17 dB. 



35 10.6 Recommended Solution 
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The recommended solution for 7 bit/s/Hz spectral efficiency is the rate 7/10 turbo code 1024 QAM 
scheme. This scheme achieves a BER of 10~ 7 at Ei/Nq = 17 dB for N=2044 information bits. 



5 11. Coding and Modulation for 12 bit/Hz Spectral Efficiency 



This section investigated one scheme that use independent I and Q modulation. The scheme 
combines a rate 12/14 coding scheme with 16384 QAM. 

10 1 LI Coding 

The coding scheme is shown in Figure 3. The two systematic recursive codes (SRC) used are 
identical and are defined in Figure 4. The co de is described by the generating polynomials 35o and 23o. 

15 11.2 Puncturing 

In accordance with an embodiment of the invention, in order to obtain a rate 12/14 code, the 
puncturing pattern used is shown in Table 22. 



20 Table 22. Puncturing and Mapping for Rate 12/14 16384 QAM 



Information bit (d) 


dl 


d 2 


d 3 


d 4 


d 5 


d6 


d 7 


dg 


d 9 


dio 


dll 


dl2 


parity bit (p) 


PI 
























parity bit (q) 














on 












128 AM symbol (I) 


(d],d2,d3, d4 d$, d6,pi) 


128 AM symbol (Q) 


(dy, dg, d9 ,dio , dn, d 12 , q7) 


16384 QAM symbol (I , Q) 


(di, d2 , d3, d4 ds, d6 , pi, dj, dg, dg ,d\0 , d\ \, d\2, qj ) 



11.3 Modulation 



25 For a 16384 QAM constellation with points at -127A, -125A, -123A, -121A, -1 19A, -1 17 A, -1 15A, -1 13A, - 

1 1 1 A, -109A, -107A, -105A, -103A, -101A, -99A, -97A, -95A, -93A, -91 A, -89A, -87A, -85A, -83A, -81 A, -79A, - 
77A, -75A, -73 A, -71 A, -69A, -67A, -65A, -63A, -61 A, -59A, -57A, -55A, -53 A, -51 A, -49A, -47A, -45A, -43A, -41 A, 
-39A, -37A, -35A, -33A, -3 1 A, -29A, -27A, -25A, -23A, -21 A, -19A, -17 A, -15A, -13A, -1 1A, -9A, -7A, -5A, -3A, -A, 
A, 3A, 5A, 7A, 9A, 1 1A, 13A, 15A, 17A, 19A, 21 A, 23A, 25A, 27A, 29A, 31A, 33A, 35A, 37A, 39A, 41 A, 43A, 45A, 

30 47A, 49A, 51A, 53A, 55A, 57A, 59A, 61 A, 63A, 65A, 67 A, 69A, 71 A, 73A, 75A S 77A, 79A, 81 A, 83A, 95A, 87A, 
89A, 91 A, 93A, 95A, 97A, 99A, 101 A, 103A, 105A, 107A, 109A, 111A, 113A, 115A, 117A, 119A, 121A, 123A, 
125A, 127 A. E av is: 
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E av =5461 A 2 (109) 

It is assumed that at time k the symbol = u 2 *S U4K U5^, u^-, uyk u$\ u ^k u\q\ U12K 
ui4^) is sent though the channel. It is assumed that at time k the symbol u\^, X12K u 3 ^ 9 u^, 115k, ug^ and u 7 ^ 
modulates the I component and ug^ 9 u ^k ? u ^k U j |k s U12K u 13k and modulates the Q component of a 16384 
QAM scheme. 

For a rate 12/14 code and 16384 QAM, the noise variance is: 



ail = 2wf^| =5461^ 



\ No 



1 ' 2 x6xEb 



No J 



No 



45 In order to study the performance of this scheme, a rate 6/7 turbo code and a 128AM is used. The 16384 

QAM scheme will achieve a similar performance in terms of bit error rate (BER) at twice the spectral efficiency, 
assuming an ideal demodulator. The puncturing and mapping scheme shown in Table 22 is for 12 consecutive 
information bits that are coded into 14 encoded bits, therefore, one 16384 QAM symbol. The turbo encoder is a rate 
12/14 turbo code, which in conjunction with 16384 QAM, gives a spectral efficiency of 12 bits/s/Hz. 

50 



45 



11.4 Bit Probabilities 
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20 



25 
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The 128AM symbol is defined as u^= (ujK i^k u$\ 114K 115k ^7^ ), where u\k is the most significant 
bit and 117k is the least significant bit . The following set can be defined. 

L bit-l-is-1 - {A64,A65>A66 , A67 , A68 , H9 » A 70 > A 71 > A72 ,A73 , A74 , A75 , A76 , A77 , A78 9 A^ Ago 
,A81 , Ag2 , Ag3, Ag4 , Agg , Ag6 , Ag7 , Agg A89 > A90 , A91 , A92 , A93 , A94 , A9g, A% ,Agj , A98 , A99 
> Aioo,Aioi,Ai02 ; Ai03>Ai04 ,A l05>Al06>Ai07,Ai08> A109, Ano .Am, 
All2,Al 13^114^115^116^117^118^119^120^121^122^123^124^125^126^127 } 

2. bit-2-is-l = {A32 ,A33 , A34 , A35 , A36 , A37 , A38 , A39 , A40 ,A4i , A42 , A43 , A44 , A45 , A46 , A47 ; A48, 
A49 , A50 , A51 , A52 , A53 , A54 , Agg , A56 ,A57 , A58 , A59 , A60 , Aei , A62 , A63, A64, A6g, A66 , A$7 , 
A68 > A 69 > A70 9 A71 , A72 ^73 , A74 , A75 , A76 , A77 , A78 , A79, Ago A81 > Ag2 , Ag3, Ag4 , Agg , A86 , 
A87 > A 88 ,A89 > A90 , A91 , A92 , A93 , A94 , A95 } 

3. bit-3-is-l = {A\ 6 ,Ai7 , Ajg , A19, A20 , A21 > A22 , A23 , A24 ,A25 , A26 , A27 , A28 , A29 , A30 , A3i ( A32 
A33 , A34 , A35 , A36 , A37 , A38 , A39 , A40 AH , A42 , A43 , A44 , A45 , A46 , A47 3 Ago ,Agi , Ag2 , Ag3, 
A84 , A85 > H6 3 Ag7 , Agg ,Ag9 , A90 , A91 , A92 , A93 , A94 , A9g ? A96 ,A97 > A98 > A99 , 
Al00 5 Al01 5 Ai02 ? Aio3 ? Al04 >Aiog,Al06>Ai07>Ai08> A 109^ Ano , Ai n } 

4. bit-4-is-l = {A8 ,A9 > A io , An , A12 , A13 , A14 , A15 , Ai<5 ,Ai7 , Aig , A19, A20 , A21 , A22 , A 2 3 , A40 
,A4i , A42 , A43 , A44 , A45 , A46 , A47 ,A4g,A49 > A50 , A51 , A52 , A53 , A54 , A55 , A72 ,A73 , A74 , A75 , 
A76 * A77 > A78 , A79 ,Ag0 A81 , Ag2 9 Ag3, Ag4 , Ag5 , Ag6 , Ag7 , A104 >Ai05,Aio6,Ai07>Ai08> Ai09, 

Ai 10 , Aj 1 1, Ai 12A1 13 M UM 1 5,Ai i6,Ai i7, A i 1&A1 19, } 

5. bit-5-is-l = {A4 , A5 , A6 , A7 , Ag ^9 9 Aio , A] 1 A20 , &2\ > A 22 > A 23 > A 24 ^25 9 A 26 , A27, A36 , A37 , 
A38 » A39 , A40 ,A4i , A42 , A43 , A52 , A53 , A54 , A55 , A56 ,A57 , A58 , A59 ; A$8 , A69 , A70 , A71 , A72 
? A73 , A74 , A7S Ag 4 , Ag5 , Ag6 , Ag7 , Ags A89 > A90 , A9i, Aioo, Aiqi, AiQ2» Ai03, A104 
,Al05,AlO6,AlO7Al i6»Ai 17^118^1 19^120^121^122^123. } 

6. bit-6-is-l = {Aq ,A\ A 6 > M > Ag ^9, A14 , Ai5,Al6 ,Al7 A22 , A23 , A24 ,A25= A30 , A3i ? A32 ,A33 A38 , 
A39 > A40 ,A4i A46 , A47, A4g,A49 A54 , A55 , A56 ? A57 A 62 » A63, A64, A65, A70 , A71 , A72 ,Aj3 > A 7S , 
A79 Ago > A 81 , Ag6 , Ag7 , Agg ,Ag9 , A94 , A 95? A96 ,A 9 7 ? Ai 0 2 ? A 10 3, A104 , Aios^ Ano , Am 

Al 12»Al 13. M\SM 19^120^121^126 > A127 } 

7. bit-7-is-l = {Ai , A2 5 A5 , A6,A9 » A io, A13 , A14A17 , Aig, A21 A22» A 25 • A26, A29 , A3o, A33 , A34, A37 
, A38A41 , A42, A45 , A46,A49 , Ago, A53 , Ag4,A57 , Agg, A^i , A62, A$5, A66, A^9 , A70A73 , A74, A77 , 
A7g ? A81 . Ag2, Agg , Ag6 ? Ag9 , A90, A93 , A94 A97, Agg ? AiOb Ai02* Ai05, A106, A 109, 
A li0» A ii3,Aii4 9 Aii7 9 Aii8,Ai21 9 Ai22 9 Ai25- A i26> 

From each received symbol, Rk , the bit probabilities are computed as follows: 



LLR( u k J = \og 



2 eXp |- 2 -' 



1__ 



R k -Ar 



S exp 

^Ajtbit-n-i s-0 



(111) 
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1 1.5 Simulation Results 
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Figure 57 shows the simulation results for 31200 information bits. A BER of 10"? can be achieved after 8 
iterations at Ei/Nq - 28.25 dB, 

12. Power vs. Bandwidth in an AWGN Channel 
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This section gives an estimate of the trade off which can be achieved between minimum required E^Nq and 
bandwidth efficiency. An information data rate of 2,044 Mbit/s and a maximum transmitter delay of 1 ms is considered. 
The corresponding interleaver size is 2,044 bits. 

12.1 Channel Model 
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All the simulations assumed the additive white Gaussian noise (AWGN) channel model, with independent I 
and Q signal. A block diagram of the system is shown in Figure 57. 

5 12.2 Simulation Results 

Simulations were run for bandwidth efficiencies from 1 to 7 bit/symbol using the recommended coding and 
modulation schemes. The results are shown in Figures 58 to 64. 

10 12.3 Conclusions 

Table 23 summarizes the minimum E^/Nq required to achieve a BER of 1 0~7. 

Table 23 Minimum Ej/Nq required to achieve a BER of 10~ 7 

15 



Spectral efficiency 


Coding Rate 


Symbol Rate 


Eb/No 




and 




For BER =10-7 


[bits/s/Hz] 


Modulation 


[ksym/s] 


[dB] 


1 


2/4 and 4 QAM 


2048 


2.2 


2 


2/4 and 16 QAM 


1024 


4.2 


3 


3/4 and 16 QAM 


682 


6.5 


4 


4/6 and 64 QAM 


512 


9.1 


5 


5/8 and 256 QAM 


408 


12.3 


6 


6/8 and 256 QAM 


342 


14.5 


7 


7/10 and 1024 QAM 


292 


17.0 



The results show the potential reduction in bandwidth for a given signal-to-noise ration for a 
particular channel. If an Eft/No of 17 dB or more is available, the symbol rate can be reduced from 2048 
ksym/s to 292 ksym/s. 

20 

13. Channel Simulation 

In this paragraph we show numerical results of the simulations for the AWGN and Impulsive nose 
using or not and outer Reed-Solomon encoder. 

25 

13.1 Without outer Reed-Solomon Encoder 
13.1.1 Net Coding Gain 

30 The net coding gain of the spectral efficiency of 4 bit/tone and 12 bits per tone, protecting more the 

information bits than the parity bits is as shown in Table 24. The case that the parity bits are more protected 
than the information bits is shown in table 25. for a BER greater than 1 0 " s the net coding gain is 0.5 dB 
worse. For BER below 10 ^ the coding gain improves respect the case that the information bits are more 
protected. 

35 

Table 24. Net Coding Gain Information bits more protected 



Bit/Tone 


Tones 


Interleaver 
Size 


# of DMT 
symbols 


Latency 
(Tx+Rx) 
ms 

< 


10- 3 


10-7 


io-9 

extrap. 




100 


5,200 


13 


10.0 


4.60 


7.42 


7.94 




800 


2 


1.5 


3,70 


4.92 


4,84 


4 




400 


1 


0.7 


3.30 


3.62 


3.84 




200 


10,400 


13 


10.0 


4.60 


7.52 


8.14 




1,600 


2 


1.5 


4.10 


6,42 


6.64 
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800 


1 


0.7 


3.70 


4.92 


4.84 


12 


100 


15,600 


13 


10.0 


4.10 


5.91 


6.03 


2,400 


2 


1.5 


3.60 


5.51 


5.63 


1,200 


1 


0.7 


3.00 


3.91 


4.03 


200 


31,200 


13 


10.Q 


4.10 


6.81 


7.53 


4,800 


2 


1.5 


3.60 


5.91 


6.43 


2,400 


1 


0.7 


3.60 


5.51 


5.63 


Table 25. Net Coding Gain parity more protected 


Bit/Tone 


Tones 


Interleaver 
Size 


# of DMT 
symbols 


Latency 
(Tx+Rx) 
ms 

< 


10- 3 


10-7 


10-9 
extrap. 


4 


100 


5,200 


13 


10.0 


4.10 


7.22 


8.04 


800 


2 


1.5 


3.20 


4.72 


5.04 


400 


1 


0.7 


2.80 


3.42 


4.04 


200 


10,400 


13 


10.0 


4.10 


7,32 


8.24 


1,600 


2 


1.5 


3.60 


6.22 


6.84 


800 


1 


0.7 


3.20 


4.72 


5.04 


12 


100 


15,600 


13 


10.0 


3.60 


5.71 


6.23 


2,400 


2 


1.5 


3.10 


5,31 


5.83 


1,200 


1 


0.7 


2.50 


3.71 


4.23 


200 


31,200 


13 


10.0 


3.60 


6.61 


7.73 


4,800 


2 


1.5 


3.10 


5.71 


6.63 


2,400 


1 


0.7 


3.10 


5.31 


5.83 



5 

13.1.2 Errors due to Impulse noise (IN) 

The impulse noise is defined as 2 consecutive DMT symbols with an increase A WGN respect to 
10 the reference noise level 10 " 7 of a carrier-to-noise ratio of 21.5 dB (spectral efficiency of 4 bit/tone) and 
45.5 dB (spectral efficiency of 12 bit/tone). 

Table 26 shows the numerical results in the case that the information bits are more protected than 
the parity bits. Table 27 shows the numerical results in the case that the parity bits are more protected than 
1 5 the information bits. 



Table 26. Error due to Impulse Noise. Information bits more protected 



Bit/ 
Tone 


Tones 


Interleaver 
Size 


RL + 
2.5 dB 


RL + 
5dB 


RL + 
7.5 dB 


RL + 
10 dB 


RL + 
12.5 dB 


RL + 
15 dB 


RL + 
17.5 dB 


RL + 
20 dB 


4 


100 


5,200 


0 


0 


0 




0 


Q 


0 


4 


800 


0 


0 


39 


65 


104 


140 


188 


243 


400 


0 


0 


10 


50 


89 


127 


161 


214 


200 


10,400 


0 


0 


0 


0 


0 


0 


0 


7 


1,600 


0 


0 


0 


127 


189 


267 


363 


448 


800 


0 


0 


40 


116 


187 


252 


346 


440 




100 


15,600 


0 


Q 


0 


0 


10 


58 


130 


207 


2,400 


0 


0 


40 


78 


121 


171 


216 


295 
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200 



1,200 
31,200 
4,800 
2,400 



43 



_75_ 
80 



98 



177 
166 



129 
90 
254 
244 



188 
175 
341 
345 



255 
313 
462 
457 



329 
482 
608 
598 



Table 27. Error due to Impulse Noise. Parity bits more protected 




13.1.3 Error Statistics 



13 1.3.1 For AWGN 

10 Table 28 shows the statistics of the errors for the case of protecting more the information bits, for a 

AWGN channel. 

Table 28. Error Statistics for AWGN. Information bits more protected 



Bit/Tone 



Tones 



100 



200 



12 



100 



200 



Interleaver 
Size 



15 



13.1.3.2 Impulse Noise 



5,20 



80 



40 



10,40 



1,60 



80 



15,60 



2,40 



1,20 



31,20 



4,80 



2,40 



1 consec. 
error 



87.30% 



94.35% 



90.28% 



89.90% 



97.94% 



90.28% 



99.79% 



98.72% 



97.94% 



99.86% 



100.00% 



98.72% 



2 consec 
errors 



10.81% 



5.64% 



9.72% 



8.63% 



3 consec 
errors 



1.47% 



4 consec 
errors 



0.29% 



0.00% 



0.01% 



1.21% 



2.06% 



9.72% 



0.21% 



1.28% 



2.06% 



0.14% 



0.00% 



1.28% 



0.00% 



0.01% 



0.00% 



0.00% 



0.00% 



0.20°/ 



0.00% 



0.00% 



5 consec 
errors 



0.03% 



0.00% 



0.00% 



0.06% 



6 consec 
errors 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 
0.00% 



0.00% 



0.10% 



0.00% 



0.00% 



0.00*3 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 



0.00% 
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Table 29 shows the statistics of the errors for the case of protecting more the information bits for a 
Impulse noise channel. 

Table 29. Error Statistics for Impulse noise. Information bits more protected 

5 



Bit/Tone 


Tones 


Interleaver 
Size 


1 consec. 
error 


2 consec 
errors 


3 consec 
errors 


4 consec 
errors 


5 consec 
errors 


6 consec 
errors 


4 


100 


5,20 


100.00% 


0.00% 


0.00% 


0.00% 


0.00% 


0.00% 


80 


75.97% 


18.99% 


18.99% 


3.36% 


0.84% 


0.84% 


40 


79.89% 


17.24% 


1.72% 


1.15% 


0.00% 


0.00% 


200 


10,40 


100.00% 


0.00% 


0.00% 


0.00% 


0.00% 


0.00% 


1,60 


80.24% 


13.05% 


4.47% 


1.68% 


0.47% 


0.00% 


80 


79.03% 


17.50% 


2.46% 


0.46% 


0.46% 


0.00% 


12 


100 


15,60 


95.19% 


4.81% 


0.00% 


0.00% 


0.00% 


0.00% 


2,40 


94.61% 


5.28% 


0.11% 


0.00% 


0.00% 


0.00°/^ 


1,20 


93.63% 


5.95% 


93.63% 


0.00% 


0.00% 


0.00% 


200 


31,20 


93.25% 


6.65% 


0.00% 


0.00% 


0.00% 


0.00% 


4,80 


94.89% 


4.95% 


0.16% 


0.00% 


0.00% 


0.00% 


2,40 


94.59% 


5.36% 


0.06% 


0.00% 


0.00% 


0.00% 



It is interesting that for the large turbo decoders the impulse errors still tends to stay within the 2 
DMT symbols. This implies a moderately large turbo coder of 5 ms follow by a convolutional 
interleaver/Reed Solomon of 10 ms should create both robust performance and good impulse resistance. 

10 

13.2 With Reed-Solomon 
13.2.1 Net Coding Gain 

15 The coding gain of the spectral efficiency of 4 bit/tone and 12 bits per tone, protecting more the 

information bits than the parity bits is as shown in Table 30. The net coding gain is shown in table 31. 



Table 30. Coding Gain with Reed-Solomon and information bits more protected 



BitTTone 


Tones 


Interleaver 
Size 


# of DMT 
symbols 


Latency 
(Tx+Rx) 
ms 

< 


10- 3 


10- 7 


10-9 
extrap. 




100 


5,200 


13 


10.0 


5.00 


8.62 


9.64 




800 


2 


1.5 


3.50 


7.12 


8.44 


4 




400 


1 


0.7 


3.50 


6.42 


7.44 




200 


10,400 


13 


10.0 


5.30 


8.82 


9.84 




1,600 


2 


1.5 


4.60 


7.72 


8.74 






800 


1 


0.7 


3.50 


7.12 


8.44 




100 


15,600 


13 


10.0 


4.40 


7.71 


8.53 




2,400 


2 


1.5 


4.60 


7.41 


8.33 


12 




1,200 


1 


0.7 


4.10 


6.81 


7.63 




200 


31,200 


13 


10.0 


4.40 


7.71 


8.53 




4,800 


2 


1.5 


4.40 


7.21 


8.13 






2,400 


1 


0.7 


4.60 


7.41 


8.33 



20 



50 



Table 3 1 . Coding Gain with Reed-Solomon and information bits more protected 



Bit/Tone 


Tones 


Interleaver 
Size 


# of DMT 


Latency 
fTx+Rxl 
ms 

< 


10 - 3 


10- 7 


10-9 
extrap. 


4 


100 


5,200 


13 


10.0 


3.42 


7.04 


8.06 


800 


2 


1.5 


1.78 


5.40 


6.72 


400 


1 


0.7 


1.94 


4.86 


5.88 


200 


10,400 


13 


10.0 


3.72 


7.24 


8.26 


1,60(1 


2 


1.5 


2.88 


6.00 


7.02 


800 


1 


0.7 


1.94 


5.56 


6.88 


12 


100 


15,600 


13 


10.0 


0.20 


3.51 


4.33 


2,40q 


2 


1.5 


0.02 


2.83 


3.75 


1,200 


1 


0.7 


-0.24 


2.47 


3.29 


200 


31,200 


33 


10.0 


1.06 


4.37 


5.19 


4,800 


2 


1.5 


0.76 


3.57 


4.49 


2,400 


1 


0.7 


0.26 


3.07 


3.99 



13.2.2 Errors due to Impulse noise 

5 

The impulse noise is defined as 2 consecutive DMT symbols with an increase AWGN respect to 
the reference noise level 10 " 7 of a carrier-to-noise ratio of 21.5 dB (spectral efficiency of 4 bit/tone) and 
45.5 dB (spectral efficiency of 12 bit/tone). 

10 Table 32 shows the numerical results in the case that the information bits are more protected than 

the parity. Table 33 shows the numerical results in the case that the parity bits are more protected than the 
information bits. 

Table 32. Error due to Impulse Noise. Information more protected 

15 



Bit/ 
Tone 


Tones 


Interleaver 
Size 


RL + 
2.5 dB 


RL + 
5dB 


RL + 
7.5 dB 


RL + 
10 dB 


RL + 
12.5 dB 


RL + 
15 dB 


RL-f 
17.5 dB 


RL + 
20 dB 


4 


100 


5,200 


0 


0 


0 


0 


0 


0 


0 


0 


80Q 


0 


0 


0 


0 


0 


0 


0 


0 


400 


0 


0 


0 


Q 


0 


0 


0 


0 


200 


10,400 


0 


Q 


0 


0 


0 


0 


0 


0 


1,600 


0 


0 


0 


0 


0 


d 


0 


0 


800 


0 


0 


0 


0 


0 


0 


0 


0 


12 


100 


15,600 


0 


0 


0 


0 


10 


58 


130 


207 


2,400 


0 


0 




0 


0 


0 


0 


0 


1,200 


0 


0 


0 


0 


0 


0 


0 


0 


200 


31,200 


0 


0 


0\ 


0 


90 


175 


313 


482 


4,800 


0 


0 


0 


0 


0 


10 


11 


65 


2,400 


0 


0 


Q 


0 


9 


10 


24 


115 


Table 33, Error due to Impulse Noise. Parity more protected 


Bit/ 
Tone 


Tones 


Interleaver 
Size 


RL + 
2.5 dB 


RL + 
5dB 


RL + 
7.5 dB 


RL-f 
10 dB 


RL-f 
12.5 dB 


RL + 
15 dB 


RL + 
17.5 dB 


RL + 
20 dB 






5,200 


0 


0 


0| 


0 


0 


0 


0 


0 



100 







800 


0 


0 


fl 

0 


ft 


ft 


A 

u 


fi 

v 


fi 


400 


0 


0 


0 


0 


0 


0 


0 


0 


200 


10,400 


0 


0 


0 


0 


0 


0 


0 


0 


1,600 


0 


0 


0 


() 


0 


0 


0 


0 


800 


0 


0 


0 


0 


0 


Q 


0 


0 


12 


100 


15,600 


V 


A 

u 


A 

u 


A 

u 


A 


A 


A 


A 


2,400 


0 


0 


0 


0 




0 


0 


0 


1,200 


0 


0 


0 




0 


0 


0 


0 


200 


31,200 


0 


0 


0 


0 


o 


0 


6 


0 


4,800 


0 


0 


0 


0 


0 


0 


0 


0 


2,400 


0 


0 


0 


0 




0 


0 


0 



13.2.3 Error Statistics. 



The statistics results obtained are practically the same than for the non Reed-Solomon case 

5 

14 Simulation using an analytical interleaver 

In this section we present the same simulation of the previous section, but using an analytical 
interleaver. 

10 

The analytical interleaver proposed here is a helical odd-even smile interleaver. 

This section gives an estimate of the trade off which can be achieved between minimum required 
E]y/No and bandwidth efficiency. An information data rate of 2, 100 kbps and a maximum transmitter delay 
15 of 1 ms is considered. The corresponding interleaver size is 2,100 bits (multiple of 2, 4, 6, 10, 14). 

14.1 Channel model 

All the simulations assumed the additive white Gaussian noise (AWGN) channel model, with 
20 independent I and Q signals. A block diagram of the system is shown in Figure 57. 

14.2 Simulation Results 

Simulations were run for bandwidth efficiencies from 1 to 7 bit/symbol using the recommended 
25 coding and modulation schemes. The results are shown in Figures 65 to 71 . 

14.3 Conclusions 

Table 24 summarizes the minimum E^/Nq required to achieve a BER of 10~ 7 . 

30 

Table 24 Minimum Ei/Nq required to achieve a BER of 10" 7 



Spectral efficiency 


Coding Rate 


Symbol Rate 


Eb/No 


n 


and 




For BER =10-7 


[bits/s/Hz] 


Modulation 


[ksym/sj 


[dB] 


1 


2/4 and 4 QAM 


2048 


2.7 


2 


2/4 and 16 QAM 


1024 


4.7 


3 


3/4 and 16 QAM 


682 


6.75 


4 


4/6 and 64 QAM 


512 


9.6 


5 


5/8 and 256 QAM 


408 


12.8 


6 


6/8 and 256 QAM 


342 


15.0 


7 


7/10 and 1024 QAM 


292 


17.5 


12 


12/14 and 16384 QAM 


170 


28.5 
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The results show the potential reduction in bandwidth for a given signal-to-noise ration for a 
particular channel If an Eft/No of 17 dB or more is available, the symbol rate can be reduced from 2048 
ksym/s to 292 ksym/s. 

5 

Due to the savings in memory and the low degradation in this kind of interleaver for a BER of 10* 
, below 0.5 dB } we recommend the use of this helical smile odd-even interleaver for these kind of 
applications. 

10 15. Conclusions 

Table 25 summarizes all simulation results from this study. 

Table 25. Summary of Simulation Results 

15 



Spectral 
efficiency 

n 


Coding 

rvuit; 


• 

Modulation 


Interleaver size 

Tri'Fni*fn<i'f"i/*in 
UUurilldvIUfl 

bits 


Required 
Eb/No [dBJ 

&-type 
interleaver 


Required 
Eb/NQ [dB] 
Analytical 
interleaver 




7/£ 
Z/O 




1 noA 


1 7 


l . / 


1 
J 


7 //I 








7 A 


L 


A if, 


o VjYVlvl 




J.J 


O. I 


7 


1/"? 


1 A A A Ayf 

ID V^/ViVl 


7^£ 


u.o 


7 7 


Z 


1 /7 




777 

LI L 


7 ft 


7 ^ 


7 

z. 


1/? 


16 DAM 


S1? 

J 




J.O 


2 


1/2 


16 QAM 


768 


4.9 


5.3 


2 


1/2 


16 QAM 


1,024 


4.5 


5.0 


2 


1/2 


16 QAM 


2,048 


4.2 


4.7 


2 


1/2 


16 QAM 


32,728 


2.9 


3.3 


3 


3/4 


16 QAM 


2,048 


6.5 


7.0 


3 


3/4 


16 QAM 


6,144 


5.6 


6.1 


3 


3/5 


32 QAM 


6,144 


6.0 


6.5 


3 


3/6 


64 QAM 


6,144 


6.1 


6.6 


4 


4/6 


64 QAM 


2,048 


9.1 


9.6 


4 


4/6 


64 QAM-1 


4,096 


8.3 


8.8 


4 


4/6 


64 QAM-2 


4,096 


10.5 


11.0 


4 


4/6 


64 QAM-3 


4,096 


1L5 


12.0 


4 


4/6 


64 QAM-4 


4,096 


11.5 


12.0 


5 


5/6 


64 QAM 


5,120 


13.0 


13.5 


5 


5/7 


128 QAM 


5,120 


13.0 


13.5 


5 


5/8 


256 QAM 


2,048 


12.3 


12.8 


5 


5/8 


256 QAM 


5,120 


11.8 


12.3 


6 


6/8 


256 QAM 


2,048 


14.5 


15.0 


6 


6/8 


256 QAM 


6,144 


14.2 


14.7 


6 


6/9 


512 QAM 


6,144 


14.5 


15.0 


7 


7/10 


1024 QAM 


2,044 


17.0 




12 


12/14 


16384 QAM 


31,200 


28.0 


28.5 



Source code for the simulations is provided in the attached Computer Program Listing Appendix. 



16. Computational Complexity saving 

To fix ideas lets use the 16 QAM case. 

Using independent I and Q with the Gray mapping, the constellation will looks like the one shown 
in Figure 13. 
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The conventional technique for extraction LLR soft-decision information from the channel is create 
a value representing the probability of the received symbol being a one as: 



X of the measures with the transmit symbol was 0 (112) 
S of the measures with the transmit symbol was 1 

where the measure is defined as: 

(~n* metric) 

io e (ii3) 

where: 

metric = Euclidian distance (or square of the Euclidian distance) from the possible transmit symbol 
15 to the received symbol. 

In Figure Receiving the point x if we do not use an independent I and Q technique, we have to 
compute the 8 exponentials (4 in each dimension) that is the distance of the point x to the four point in the I 
dimension and to the 4 points in the Q dimension as shown (xl, x2, x3 5 x4, yl, y2, y3, y4). To compute the 
20 4 LLR probabilities we need 4* 14 adds and 4* 16 multiplications and 4 divisions. The total number of 
operations is 132. 



25 
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With the independent I and Q if we want to compute the LLR of the least significant bit, 



X Jxhy]) + e (xj+y 2 2 ) + e (xl+y 2 3 ) + ^fry 2 4 ) + Jfry\) + ^x\+y\) + ^fry\) + e (xW 4 ) 



bit=0 



(114) 



£ e (*W<) + e (*W 2 ) + e (xW 3 ) + e (xl+y 2 4 ) + Jxfry 2 ,) + JxWz) + e (xW 3 ) + e (xW 4 ) 

(e*' + e x ')e (/ ' +/ > +y > +y * J _ e x '+e xi 
( e x\ + g xi ) e (yWfy>y\) gd + e * 2 ' 

With the independent I and Q if we want to compute the LLR of the second-least significant bit, 

Z 

X (e 4 + e x ')e (/ ' +y ' +y ^ } e 4 + e x 



bit=0 



With the independent I and Q if we want to compute the LLR of the third-least significant bit, 

35 



Z r^+e*v**** (116) 

bit=0 
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With the independent I and Q if we want to compute the LLR of the most significant bit, 



£ (e y2 ^^)e ( ^ x¥x¥A) e yl +e y < 



(117) 



bu=0 



The same reduction of calculations occurs for all bits. 

The number of exponentials to compute is the same, 8, the number of additions is 4*2=8, the 
10 number of multiplications is 0 and the number of divisions is 4. The total number of operations is 20. 

When the constellation size increase the saving in computations also increase considerably. 

For the case of 16 QAM, 64 QAM, 256 QAM and 16384 QAM the comparison is as follows: 

15 



16 QAM 


#exp 


#add 




#div 


TOTAL 


Full 


8 


4*14 


16 


4 


84 


Independent I&Q 


8 


4*2 


0 


4 


20 



20 



64 QAM 


#exp 


#add 


#mul 


#div 


TOTAL 


Full 


16 


6*62 


64 


6 


458 


Independent I&Q 


16 


6*6 


0 


6 


58 




128 QAM 


#exp 


adds 


mul 


div 


TOTAL 


Full 


32 


7*126 


121 


7 


1,042 


Independent I&Q 


32 


7*14+2*14 


14 


14 


186 



256 QAM 


#exp 


#add 


#mul 


#div 


TOTAL 


Full 


32 


8*254 


256 


8 


2,328 


Independent I&Q 


32 


8*14 


0 


8 


152 



16384 QAM 


#exp 


#add 


#mul 


#div 


TOTAL 


Full 


256 


14*16382 


16384 


14 


246,002 


Independent I&Q 


256 


14*126 


0 


14 


2,034 



25 

Nodd 



N=2 n QAM 


#exp. 


#adds 


#mul 


#div 


TOTAL 


Full 


2(2N)"2 


n(N-2) 


N-n 


n 


2(2N)l'2+(n+l)N- 
2n 


Independent I&Q 


2(2N)l/2 


2n 2 +4n 


2n 


2n 


2(2N)i/2 + 2n2+8n 



30 



55 



N even 



N=2 n QAM 


#exp 


#add 


#mul 


#div 


TOTAL 


Full 


2N1/2 


n(N-2) 


N 


n 


2Nl/2 + ( n +i)N-n 


Independent I&Q 


2N"^ 


n(Nl^-2) 


0 


n 


(n+2)N 1 ^+n 



5 The increase complexity for this type of constellation can be shown to be of 0((N) i/2 ) where N is 

the number of constellation points. For the odd constellation the procedure defined in Provisional Patent 
Application Serial No. 60/248,099 has been used. 



17. Use of Puncturing Patterns and Protection on Parity/Information Bits 

10 

In conventional Turbo Codes, when no comment about puncturing is made, it is understood that the 
puncturing pattern used selects one parity bit from the first encoder as one encoded bit and parity bit from 
the second encoder as the next encoded bit in an alternating fashion. Thus in the conventional Turbo Code 
using a four bit constellation, two parity bits are selected for every two information bits, and as such the 
1 5 number of parity bits represented in each symbol is equal to the number of information bits represented in 
the symbol. 

In accordance with embodiments of the invention, for the constellations selected, one or two parity 
bits are sent in each dimension, providing, for this reason, the possibility to protect the parity bit(s) more 
20 than the information bit(s), or less than the information bit(s), because of the Gray mapping. The number of 
parity bits represented in each symbol is less than the number of information bits represented in the symbol. 

In an AWGN channel, if the information bits are more protected, for BER higher than 10" 7 , the 
performance is 0.5 dB better. This means that for the same BER, the EJN^ in the case that the information 
25 bits are more protected, is 0.5 dB lower than the EJN 0 if the parity bits are more protected. For BER below 
10" 7 , this statement is reversed. For the same BER, the EJN 0 , in the case that the information bits are more 
protected, is 0.5 dB higher that the Et/No if the parity bits are more protected. 

In an Impulse Noise environment, the greater protection of the parity bits provides more 
30 immunization against impulses (around 5 to 10 dB more). Once the impulse affects the more protected 
parity, they are more harmful than the case where the parity is less protected and the information is more 
protected. 

Figure 74 shows the performance of the 6/4 rate 64 QAM for the case that the parity bits are most 
35 protected and for the case that the parity bits are least protected. 

A detail explanation of the encoding and puncturing procedures is shown in Figure 75, where 
matrix |A|, is the total set of bits that are encoded. The number of rows of this matrix is the number of input 
symbols to be transmitted, "m". The number of rows of this matrix is the number of bits encoded in each 
40 symbol, "n'\ 



\A\ = 



m,n j 
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Each row of the matrix |A| is sent to the Turbo Encoder once every symbol time. For symbol j, the 
input to the turbo encoder is |A j| = {a - h x , a j2 a j>n }. For each symbol, the output of the first 
convolutional encoder will be |P j| = p jt 1 , p h2 p j>n , and the output of the second convolutional encoder 
will be |Q j| = q j,i ? q j, 2 q j, n . Each encoded bit a produces one parity bit in the first encoder p M , 
5 and one parity bit in the second encoder, q Essentially the encoder runs w n" times for each symbol. 



PU P\,2 



Pu 



\p\ = 



K Pm,\ Pm,2 • 



#1,2 



; Pj = {Pj,x,Pj,2, ,Pj,J o»9) 



PmjiJ 
$0 



101 = 



; Q J ={<i j ,x,q h i, ,q } ,} (120) 



10 



15 



<^7m,1 <7m,2' 



1m,n , 



The output of the first convolutional encoder and the output of the second convolutional encoders 
go to the puncturing pattern block, as shown in Figure 75. This puncturing pattern block, selects the parity 
bits to be transmitted. The first output provides the parity bits to be transmitted from the first encoder |P' j| = 
P' i, i » P' j,2 >— » p' j, a • Th e second output provides the parity bits to be transmitted from the second encoder 

IQ'jl = q'j,i>q'i,2v..,q' j , (J . 



P\,l Pl,2 ■ ■ P\ A 



\P m,] P m,2 



Pm,a 



; Fj = {p , JA ,P , j*, ,p'jJ (i2i) 



\Q\ = 



9 m,1 Q m,2 



20 



; Qj = W 3 ,x,<i , j ,2> 022) 



In accordance with an embodiment of the invention, a + P < n. The total number of parity bits sent 
in each symbol is a + (3, In traditional puncturing, a + (3 = n and a = (3 — n/2. 



The implementation shown in Figure 76 includes additional uncoded information bits, where 
25 matrix |B| is the total number of information bits. This is the general form for Multi-level Turbo Codes, 
where not all the information bits are coded. The number of rows of this matrix is the number of input 
symbols to be transmitted, "m". The number of rows of this matrix is the number of bits in each symbol, 
a N w . Matrix |B| has a sub-matrix |A|, which represents the information bits to be encoded. The number of 
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rows of this matrix is the number of input symbols to be transmitted, "m". The number of rows of this 
matrix is the number of bits that are encoded in each symbol, "n'\ 



10 



B\ = 



1JV 



18. Puncturing Rate Adaptation 



a, 



a. 



1,2 



a, 



In 



J m,N J 



(123) 



Using the Channel State Information (values of gi and bi for each tone for the case of a G.992.1 or 
G.992.2 modem), the puncturing and mapping can be adapted to alter the transmitted signal bandwidth, 
which allows communication at a desired BER. 



If there are too many errors (i.e., BER is too high), this can be improved by increasing the 
constellation size and reducing the puncturing, thereby sending more parity bits. The puncturing of the 
system is reduced, increasing the values of a and p. If the values of a and p are increased to twice their 
15 original values, the redundancy of the communications is also increased to twice its original value. 

Alternatively, if there are too many errors (i.e., BER is too high), this can be improved by reducing 
the information data rate and reducing the puncturing, thereby sending more parity bits. The puncturing of 
the system is reduced, increasing the values of a and (3. If the values of a and P are increased to twice their 
20 original values, the redundancy of the communications is also increased to twice its original value. 

If the channel is working with a very low number of errors, it is possible to send information with 
fewer parity bits, by increasing the puncturing which reduces the values of a and p. 

25 In the case that the bandwidth of the system is not a concern (i.e. spread spectrum systems or low 

data rate situations), reducing the puncturing increases the performance. For this case ct and p should use 
their maximum value, "n". 
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The monitoring of the error rate may be performed in the receiver using various techniques. Noise 
measurements based on the distance between each actual QAM signal point received and its nearest 
constellation points is a frequently used MSE (Mean-Squared Error)performance metric. Alternatively, the 
detection of errors in a bit-stream by a forward error-correction protocol, such as Reed Solomon or Turbo 
Decoder, may be used to determine performance metric. Higher-level protocol detection of errors, such as 
from an HDLC CRC checker, may also be used to determine a performance metric. Thus the actual error 
rate need not be measured. Rather, a performance metric may be derived in the receiver using any suitable 
technique and then characterized to determine when the receiver would benefit from a change in the 
puncturing pattern. Once the receiver decides to select a different operating puncturing pattern, it 
communicates this selection to the transmitter using a renegotiation process. 

Taking into account these possibilities, it is clear that the system may adapt to different scenarios 
modifying the puncturing pattern using a "puncturing pattern adaptation" based on the condition of the 
channel, increasing or reducing the values of a and p. The operation of the channel can be controlled using 
a a-mask and P-mask. The a-mask and p-mask are V bits long; each position can take a "0" value or a 
"1" value. 

The "0" value in the position "i" in the a-mask means that the puncturing pattern will not include 
the parity bit "p h " from the first encoder. The "1" value in the position "i" in the a-mask means that the 
puncturing pattern will include the parity bit "p p " from the first encoder. The "0" value in the position "i" 
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in the P-mask means that the puncturing pattern will not include the parity bit "q h " from the second 
encoder. The "1" value in the position "i" in the p-mask means that the puncturing pattern will include the 
parity bit u q from the second encoder. Equations (124) and (125) shows the values of these masks. 

5 

ajinask= {a Jly a jf2 ,"; a Jn } a Jt = 0,l i = L...n (124) 
fijmask= {fi u ,fi u ,... 9 fi u } fi^O,! i = 7..../i (125) 

10 

As seen in the puncturing patterns disclosed in Tables 7-22, the a-mask and P-mask of a 
puncturing pattern typically utilize the same repetitive bit selection pattern, with the pattern in theP-mask 
being shifted by one or more bits with respect to the pattern in the a-mask. The patterns of the a-mask and 
P-mask are typically non-coincidental, in that the two patterns do not select a parity bit at the same time. 

15 

19. Summary of operation of transmitters and receivers in preferred embodiments of the invention 

Figure 77 illustrates a basic process in a transmitter in accordance with an embodiment of the 
invention. The process produces a modulated signal with forward error correction from an information bit 

20 stream in a QAM transmitter. The transmitter produces 100 parity bit streams that correspond to an inputted 
information bit stream using first and second concatenated coders interconnected by an interleaver. The 
transmitter selects 102 subsets of the first and second parity bit streams in accordance with a puncturing 
pattern. The transmitter combines 104 the selected subsets of said first and second parity bit streams with 
said information bit stream. A QAM symbol stream is produced 106 by mapping a first subset of the 

25 combined bit streams to an I dimension and mapping a second subset of the combined bit streams to a Q 
dimension. The QAM symbol stream is modulated 108 to produce a modulated signal, and the modulated 
signal is transmitted 1 10 over a communication link. 

As described above, any puncturing pattern may be implemented in accordance with the above 
30 process. However, it is preferred to implement one of the novel puncturing patterns disclosed herein. As 
described above, the novel puncturing patterns disclosed herein are characterized in that, for an c-bit 
constellation, wherein each symbol represents n information bits and c - n parity bits (also a+ p parity bits 
using the terminology of Section 17 above), the number of parity bits c-n selected by the puncturing 
pattern is less than n. In other words, the number of parity bits represented in each symbol is less than the 
35 number of information bits represented in that symbol. Examples of such puncturing patterns are provided 
in Tables 7 and 9-22, and those of ordinary skill will recognize a variety of other patterns that may be 
implemented in accordance with the parameters described herein. 

Figure 78 illustrates a basic process in a receiver that is complementary to a transmitter as 
40 described above. The process recovers an information bit stream from a noisy modulated signal with 
forward error correction in a QAM receiver. The receiver receives 120 a modulated signal from a 
communications link. The signal includes errors. The received signal is demodulated 122, and a decoded 
bit stream is produced 124 by iteratively decoding the demodulated signal using first and second 
concatenated coders connected by an interleaver. The demodulated signal is decoded independently in the I 
45 dimension and the Q dimension using a puncturing pattern. The puncturing pattern corresponds to a 
puncturing pattern implemented in the transmitter that produces the modulated signal. The receiver 
regenerates 126 the information bit stream from the decoded bit stream. 

As noted above, embodiments in accordance with the invention may adapt the puncturing pattern 
50 used in the transmitter in accordance with conditions in the communication system. Figure 79 illustrates a 
basic process in a receiver for adapting to conditions in the communication system. The receiver receives 
130 a QAM symbol stream from a transmitter of the communication system. The receiver determines 132 a 
performance metric of the communication system based on the received QAM symbol stream. The 
performance metric may be determined based on detection of errors in the QAM symbol stream in 
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accordance with the error detection techniques discussed above or any other conventional error detection 
technique. If the performance metric is unsatisfactory, the receiver selects 134 a new puncturing pattern for 
use in the transmitter to improve the performance of the communication system. The puncturing pattern is 
selected based on the performance metric. The puncturing pattern is then communicated 136 to the 
5 transmitter. This may be accomplished by communicating the actual pattern to the transmitter, or by 
transmitting an identifier of a predefined puncturing pattern. 

While the foregoing disclosure addressed a variety of preferred embodiments of the invention, 
those having ordinary skill in the art will recognize a variety of further embodiments that may be 
1 0 implemented in accordance with the novel principles disclosed herein. 
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